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Preface 



Although Charles Babbage laid down several ideas for computing "en- 
gines," the forerunners of today's computers were largely developed in the 
early 40's as part of the war effort. The Robinson series cryptoanalytic ma- 
chines developed in England in 1941 spawned many families of computers 
still in use today. The MIT differential analyzer and real-time aircraft simula- 
tion project led to the Whirlwind, and eventually to the immensely successful 
DEC family of PDP computers (Programmed Data Processors). And, of 
course, Eniac built at the University of Pennsylvania, was the guiding light 
behind Univac, IBM and many other successful manufacturers. 

Many people today poke fun at these early machines and regard them as 
dinosaur-like relics. However, it is interesting to consider that a large-scale 
computer of about 25 to 30 years ago had about the same amount of power as 
a typical personal computer of today. It was generally not as reliable or user- 
friendly as a personal computer, and, of course, cost tens of thousands times 
as much. Why bring this up? 

Because a computer of the 50's required that the programmer be very 
clever and resourceful to solve problems within the capabilities of the com- 
puter. He did not have vast gobs of memory available, blinding quick calcula- 
tion speed, or random disk access. In other words, he had about the same 
problem to face as you do with your personal computer. 

I do not mean to imply that your personal computer is not a full-fledged 
computer. It certainly is just as much a computer as a room-filling giant of to- 
day. However, because of the relatively small memory, it cannot store a large 
data base. Nor is it suitable for extensive word processing or massive calcula- 
tions. Highly detailed graphics are best left to other machines as well. 

What can we learn from the computing pioneers of the 50's that will help 
us today? Perhaps most important is the discipline of thoroughly analyzing a 
problem, breaking it down into manageable steps, and solving it a step at a 
time. It is also important to determine what can be done "off line" and what 
must be done on the computer. 

IX 



That is what this book is all about. While it has more than 50 ready-to-run 
programs, the main thing you should look for from the book is an approach 
to solving problems — big and small. Some of the problems demonstrate the 
capability of the computer; others identify its shortcomings. It is important to 
be familiar with both the strengths and weaknesses of your tools so you can 
recognize the types of jobs for which they are suitable (and not suitable). 

This book focuses primarily on mathematical and educational applications 
for the computer. There are many other excellent sources of information 
about other applications and for making best use of your personal computer. 
Using the approaches described in this book should enable you to easily con- 
vert programs and use applications from other books and magazines such as 
Creative Computing. 

This book is designed to be read with a working computer at hand. While 
there is textual material to be read, the most important things are the experi- 
ments and problems to be tried with your own computer. The book raises 
many questions for which you should try to find answers. There are no an- 
swers to these questions and problems in the back of the book; you should be 
able to discover the answers as you work the problems out on your computer. 

You will be able to incorporate many of the routines and approaches in the 
book into programs of your own as you use your computer to deal with "real 
world" problems. Other programs will simply point you in the right direc- 
tion. And some of the programs in the book are just plain fun. Learn. Experi- 
ment. Have fun! 

Morristown, New Jersey 

March 1983 David H. Ahl 
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Drill and Practice 



Throughout hfe, there are certain things that simply must be memorized. 
Obvious things that fall into this category are the addition and multiplication 
tables, the spelling and meaning of words, how to tell time, and the monetary 
system. 

However, depending upon one's chosen profession, there are many other 
things to be memorized. A doctor must know what diseases match what 
symptoms. A chemist must know the gas laws, the properties of elements and 
so on. A pilot must instantaneously know the meaning of readings on scores 
of instruments. 

To memorize a set of facts, you must go over them again and again and 
keep trying different variations. Here is where the computer comes in. It is 
able to present randomly scores of different problems to you for as long as 
you wish. Some programs will automatically adjust to your level of com- 
petence and will grade you; other programs simply present the problems and 
leave the grading up to you. 

There are four programs in this chapter and two in the Science chapter 
which present material in a drill and practice format. Examine the methods 
used in these programs and then make up some drill programs of your own 
for subjects with which you are having trouble, or make up programs for 
other members of your family. 
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Addition Practice 

This program demonstrates a simplified addition drill and practice routine. 
This type of drill is sometimes called computer assisted instruction (CAI), al- 
though CAI can also apply to tutorial and other approaches as well. 

When the program is run, it will first ask "No. of digits?" You enter a 
number and each addend will contain that many or fewer digits. 

The program will present any number of practice problems that you spec- 
ify. The program presents each problem in turn. The program will not pro- 
ceed to the next problem until the current one has been answered correctly. 
After the last problem is answered correctly, the score is printed with an 
appropriate comment. 

There are many improvements and extensions possible in a program like 
this one. For example, you might want to modify the program so it tells the 
user the correct answer after a problem has been answered incorrectly two (or 
three) times. 

A more complicated modification would be to change the program to 
present different kinds of arithmetic problems such as subtraction, multiplica- 
tion, and division. 

Some of these modifications have been made in the next program. 
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Addition Practice, 
Adjusted by Grade Level 

One of the major disadvantages with many drill and practice exercises is 
that they tend to be either boring or frustrating, depending upon the ability of 
the user relative to the level of the material. To compensate for this, a method 
is needed which will adjust the difficulty of the problems to the ability of the 
user. 

Ideally, such a system would weigh the most recent performance most 
heavily but would not ignore previous performance. It should allow a user to 
advance to more difficult problems than his current mastery level. It should 
also continue to give some practice on problems already mastered. 

Some commercial software packages approach these goals along traditional 
lines, i.e., determine in which type of problems a student should receive prac- 
tice by using a complicated computer managed instruction score recording 
and adjustment system. 

The approach here is more innovative; it uses a single measure for each 
type of problem — call it "estimated grade level" — which meets all of the 
objectives stated above. 

How does it work? The most recent problem presented counts 10% of the 
overall score if it was answered correctly and was over the current user grade 
level, or if it was answered incorrectly and was under the current user grade 
level. Otherwise it is ignored. This may be easier to visualize in the form of a 
chart: 

Answer 
Right Wrong 



Higher than Raise student 
grade level grade level 



Ignore 



Problem 



Lower than Ignore Lower student 

grade level grade level 



At first glance this might look complex and somewhat goofy, however, 
what it really means is that a student is rewarded for doing a problem beyond 
his grade level but he is not penalized if he cannot do it. On the other hand he 
is penalized if he cannot do a problem lower than his grade level, but is not 
rewarded for doing one lower. 

Each problem affects the estimated grade level a little bit, with the most 
recent problems being weighed the most heavily. If the current grade level of 
a student is 1. and the level of the most recent problem to be averaged in is P, 



then the averaging formula is simply: 

L = .9L + .IP 

The remaining task before a program can be written is to assign a grade 
level to each problem presented. Unfortunately, this will vary depending 
upon the local school system, the textbook used, and the teaching method. 
Also a huge data base can not be stored in a small computer, so it is desirable 
to devise a simple method of determining grade level for different problems. 
One straightforward approach is to present problems up to one-half a grade 
level over and under where the student currently is. Thus the overall range of 
problems for a student at grade level 3.2 would be 2.7 to 3.7. 

How do we generate the right problems? Consider one type of skill, vertical 
addition. It is normally introduced in the first grade and continues through 
Grade 4 (actually 4.9). The simplest problem in this program is 1 + 1 and the 
most difficult is 999 + 999. Since learning is not a linear process (it is slow at 
first, and then progresses rapidly), an exponential formula can be used. For 
example; 

Addend = 1.73 x (Grade level)" 
or 

Grade level = V Addend/ 1.73 

This means that students at various grade levels will be working with the 
following maximum addends: 

Grade level Addend 

1.0 1 

2.0 27 

3.0 140 

4.0 442 

4.9 997 

Now it is a relatively straightforward, although somewhat tedious, matter 
to tie all these elements together in a computer program. 

A few notes about the program. The variable G2 is the problem grade level 
that is always within one-half of a grade level of the current student level, Gl. 
The complicated mess in Statement 340 produces a moving grade level 
average. 

The recording of the grade level and carrying it over to the next lesson is a 
manual process. On a computer system with a permanent mass storage de- 
vice, this would be kept on the system. 

There are many possible changes and extensions to this program. For 
example it could present different types of problems such as horizontal addi- 
tion, vertical and horizontal subtraction, as well as multiplication, division 
and fraction problems. 
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Time/Speed/Distance Problems 

As well as being able to present simple numerical drill and practice prob- 
lems, the computer can present word problems for solution as well. 

In this program the formula relating time, speed and distance was applied 
to a problem involving both a car and train. 

The problem can be stated as follows. A car traveling C miles per hour 
(computer generates an integer 40 through 65) can make a certain trip in D 
hours (computer generates an integer 5 through 20) less than a train traveling 
at T mph (computer generates an integer 20 through 39). How long does the 
trip take by car? When the two simultaneous equations are solved they pro- 
duce the single equation for the answer shown in Line 1 10. 

Notice the calculation in Line 120. This calculates the percent difference 
between the actual answer and the one entered by the user. 

Notice also that the computer calculates the correct answer in Line 1 10 and 
prints it (on the screen) in Line 200. This answer may have many decimal 
places; as the program is written it is rounded off to two decimal places. If the 
0.5 was not added in Line 110, the number would be truncated and not 
rounded off. This is an important calculation and one which you will find in 
many other programs throughout the book. 

Consider other problems that can be used as the basis for this kind of drill 
and practice exercise. Teachers, for example, might wish to have students 
write drill and practice programs on their own. Different problems could be 
given to one or a small group of students to serve as the basis for a program. 

After the programs are written, students can try out the programs of other 
class members. This approach ensures that each student not only understands 
the type of problem assigned to him, but also gets practice in solving other 
problem types as well. This is an effective technique for stimulating interest as 
well as for learning how to solve word problems. 
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Kinematics Problems 

Kinematics relates to the dynamics of motion of bodies apart from consid- 
erations of mass and force. Like many other types of problems, these can be 
generated and presented by the computer to provide practice in solving them. 

This program presents a simple kinematics problem for solution. The com- 
puter generates a new value for each problem. The problem is as follows. A 
ball (or any other object) is thrown up at velocity V meters per second (com- 
puter generates an integer between 5 and 40). The user must then calculate 
three factors about the resulting flight of the ball: maximum height, time until 
it returns, and velocity after T seconds (computer generates a time less than 
the total flight time). 

The key benefit in using a computer to present problems of this type is 
motivation. The calculations required of the user are no different than those 
in the back of a chapter in a book or those on homework assignments. How- 
ever, answering them when they are presented by the computer seems to 
make it more of a challenge and, frankly, more fun. 

The computer program checks each of your responses to see if it is within 
15% of the correct answer. If it is, your answer is considered correct. You 
may wish to change this percentage to require more or less accurate calcula- 
tions. You may also wish to change the computer calculations to round off to 
one or two decimal places. 
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"This nano-computer is great, but working the keyboard is a real problem.' 
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Problem Solving 



In many courses in school, the textbooks present a great variety of devices 
for solving the problems that have been neatly grouped together at the end of 
each chapter. Typically these devices consist of formulae, equations, rules and 
theorems. After a careful study of these devices, teachers give exams which 
test your ability to recall them. 

But what do you do if you are faced with the more realistic situation of not 
being told what device is likely to solve which problem or, worse yet, of hav- 
ing forgotten how to use a technique altogether? Is all hope lost? Of course 
not, although some people seem to believe that it is. 

In this chapter, several of these nasty devices mentioned above are pre- 
sented. For example, there are devices written into computer programs that 
will solve a quadratic or exponential equation and others that will calculate 
the roots and draw a plot of any function. 

However, you must remember that while these devices are useful in solving 
certain problems, the really important thing is to understand the underlying 
logic and approach. Then when it comes time to solve real world problems 
you will be better prepared to face them. Incidentally, many of the methods 
and devices presented in this chapter are used in later chapters to solve other 
kinds of problems. 
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How Many Tickets? 

Here is a problem. At a school raffle to raise money, the organizers have as 
a prize an electronic game for which they paid $18.00. To add interest to the 
raffle, the organizers have decided to sell tickets for an amount (in cents) 
equal to the number on the ticket for tickets numbered 1 to 50. For ticket 
numbers over 50, the price is 25 cents each. The organizers want to know 
how many tickets they must sell to exactly break even. 

It is probably easiest to visualize a problem of this sort with a flowchart. In 
the flowchart, T will equal the total money collected and will increase as 
more tickets are sold. The ticket number is N. When T equals or exceeds 
$18.00, N will be the answer. 

Note that the flowchart has two logical branching points (if statements in 
the program). The first compares the current ticket number to 50; if it is less, 
the ticket number is added to the total whereas if it is greater than 50, the 
total is increased by 25 cents. 



T = 



T=TvM 




T=T+T.S 



Nei 



PRlNr 
M,T 





NEXT N 
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The second branch point compares the total amount collected, T, to $18.00 
(actually 1800 cents). If T is equal to or greater than 1800, the break even 
point has been reached and the values of N and T are printed (on the screen). 

A problem of this kind can be done by hand, however, because of the repet- 
itive additions it is quite tedious. Also, doing it by hand frequently leads to an 
answer of 72 rather than the correct answer of 7 1 . Try it yourself and see 
what you get. 



AS ■'■ = 

20 FOR r4=l TD 2000 

3 IF N < = 5 THEN 60 

40 T=T+25 

5 li D T D 7 

6 T = T •»- H 

70 IF 7=13 00 THEN 9 

8 NEXT r^ 

90 PRINT N; "TICKETS SOLD TO 
100 PRINT "COLLECT S'-JT-'IOO 

9 9 9 END 
> 

> 

> R U N 

71 TICKETS .9DLD TD 
COLLECT J 18 

** DONE ** 



Many problems can be solved quickly and correctly with a computer using 
logical analysis and a flowchart. More complex problems may have to be bro- 
ken down into additional steps and require a longer flowchart, but the ap- 
proach is fundamentally the same. 

Here are two problems for you to solve. 

The diameter of a long-playing record is 12 inches. The unused center has a 
diameter of 4 inches and there is a smooth outer edge 1/2 inch wide around 
the recording. If there are 9 1 grooves to the inch, how far does the needle 
move during the actual playing of the recording? 

A movie theater charges $2.50 for an adult admission and $1.00 for a child. 
At closing, the cashier counted 385 ticket stubs and had $626.50 in cash. 
How many children entered? 
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Drinking and High Blood Pressure 

This program illustrates how several simple equations can be put in a com- 
puter program to solve a more difficult overall problem. 

Here is the problem. In a survey of 1000 adults, it was found that 35 had 
high blood pressure. Of those with high blood pressure, 80% drink 1 5 oz. or 
more of alcohol per week. Of those without high blood pressure, 60% drink a 
similar amount. What percent of drinkers and non-drinkers have high blood 
pressure? 

This problem requires the solution of several simple equations. They could 
all be combined into one large equation, but it may be easier to understand 
the approach (and change variables later on) by writing a program with five 
separate equations. 

If H equals the number of people with high blood pressure, then H = 35 
(Line 10). Letting HI equal the number of people with high blood pressure 
who drink leads to HI = .8 x H (Line 20). 

Letting LI equal the number of people with low blood pressure who drink 
yields LI = .6 x (1000 - H). Then, the total number of drinkers, D = HI + 
LI. 

Finally, the percentage of drinkers with high blood pressure is X = H 1 x 
100 /D. 

The program solves the problem in a jiffy. The solution for this type of 
problem can be easily written directly in Basic without any need for a flow- 
chart or detailed analysis. Recognizing this type of problem readily will save 
a great deal of pencil pushing time. 

10 H = 3 .5 

2 H 1 = . 13 * H 

3 L 1 = . 6 * < 1 - H > 
40 D = H 1+L I 

5 X = I N T <: I * < H * 1 .' D ;> > / 1 

60 E=H*100/1D00 

70 PRINT TflB<8> ; "PEOPLE WITH 

HI 6H" 

SO PRINT T fl B < 9 ;■ ; " B L D D D PRESS 

URE " 

9 PRINT "DRINKERS " ; X ; " 
'A " 

100 PRINT "ALL PEOPLE "!B 
.999 END 

>RLIN 

PEOPLE IJITH HIGH 

BLOari PRESSURE 

DRINKERS 5 . 7 ■< 

ALL PEOPLE 3.5 

** DONE ** 
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Here is a problem that doesn't require a single equation but makes use of 
the approaches discussed so far in this chapter. Can you solve it with three 
Basic statements? 

In early January, a shopkeeper marked down some calendars from $2.00 to 
a lower price. He sold his entire stock in one day for $603.77. How many did 
he have? 

Here's another easy one. A town in India has a population of 20,000 peo- 
ple. Five percent of them are one-legged and half of the others go barefoot. 
How many sandals are worn in the town? 
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Two Simultaneous Equations 

So far in this chapter, only problems with linear equations have been 
considered. But the computer can be used to solve much more difficult equa- 
tions. In fact, it is in problems involving second and third degree equations, 
exponentials, and the like where the computer really starts to pay off. Con- 
sider the following two simultaneous equations: 



2X = ibY 



16y 3^ = 27y 



It is not at all easy to solve these two equations by hand. But a simple Basic 
program can be written to solve the equations using trial and error. This is 
sometimes referred to as a brute force approach because every possible 
combination of numbers between an upper and lower limit is tried until a 
solution is reached or until the program runs out of values. 



10 FDR X=l TO 10 

2 FOR ■■!' = i T D 100 

3 IF 2 " X < > 1 6 * V -^ 3 THEN 70 

4 IF 3 - K !. > 2 7 * N' THEN 70 

4 5 PRINT " F I N ft L L Y SOLVED IT. 

II 

5 PRINT " X = " 5 X , " V = " ; Y 
€.0 STOP 

7 NEXT V 

8 NEXT X 

9 PRINT "Na INTEGER SOLUTID 
N " 

100 PRINT "BETUIEEN 1 AND 100 



9 9 9 END 



■RUN 

FINALLY SOLVED IT... 
X= 4 V= 3 

** DONE ** 



In this particular case, a solution is reached rather quickly with x = 4 and 
y = 3. However, if in the second equation the y coefficient is changed slightly 
from 27 to 28, the computer will try 10,000 possible solutions before finally 
concluding that no integer solution exists — at least within the range of to 
100. Warning: this will run for a very long time. 
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>40 IF 3-^X<>28*Y THEN 70 

> 

> R U N 
ND INTEGER SDLUTIDN 
BETWEEN 1 FIND 100. 

** DDNE ** 



Although the trial and error (brute force) approach is widely used, it is 
highly inefficient. In general, a systematic or guided trial and error approach 
is preferable to one that simply tries every possible solution. However, for 
some problems the simple "try every value" approach may be appropriate. (A 
comprehensive discussion of trial and error approaches can be found on pp 
36-40 of Computers in Mathematics: A Sourcebook of Ideas.) 

Three problem solving approaches have been discussed so far. Remember- 
ing them, how would you do this problem? A boy and his sister visited a farm 
where they saw a pen filled with pigs and chickens. When they returned 
home, the boy observed that there were 18 animals in all, and his sister re- 
ported that she had counted a total of 50 legs. How many pigs were there in 
the pen? 









^r^ll.-ill 
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Quadratic Equation Solver 



For any values A, B, and C of a first degree quadratic equation (Ax^ + Bx 
+ C = 0), this program will compute the roots of the equation. The solution 
is based on the quadratic theorem which solves for roots with the following 
formula: 



X 



-b±-\/bmac 



2A 



Assuming that A, B, and C are real numbers, the following principles 
apply: 

1 . If B^ - 4AC is positive, then the roots are real and unequal. 

2. If B^ - 4AC equals 0, then the roots are real and equal. 

3. If B^ - 4AC is negative, then the roots are imaginary and unequal. 
The program takes into account all these possibilities and correctly identi- 
fies the type of roots along with their values for any set of coefficients. 

Is this program useful by itself? Except for solving quadratic equations for 
algebra class, probably not. However, as a routine in a larger program to 
solve quadratic equations that might be encountered, it could be very useful. 



10 

3 a 

2 
DLL 

3 
40 

5 

6 
70 
80 
90 
100 
1 10 
12 
130 
140 
150 
160 
170 

18 

19 
200 
205 
£1 
R ) .-' 
220 

R y -^ 

999 



P R I r>IT 
L 'V E R " 
F' R I N T 
D LI I N 6 : 
INPUT 
INPUT 
INPUT 
PRINT 
R = B '^ 2 - 
IF fl <:: > 
F' R I H T 
P R n>i T 

STOP 
P R I fj T 
IF R <: 
P R I \\ T 
IF R = 
PRINT 
PRINT 
STOP 
P R I rj T 
STOP 
P R I tJ T 
P R I tJ T 
2 * R ; " * 
PRINT 

2 # R ; " # 

END 



UftD 

LEft 

ALU 
ALU 
RLU 

R*C 
THE 
HRT 
DEG 

THE 
THE 

THE 
-E + 
-E- 



R R T I C E Q U fl T I a N 
SE ENTER THE F 



E DF R=" 
E D F B = " 
E D F C = " 



N 1 2 

IS R FIRST" 

REE ECURTIQN" 

R D a T S R R E " ; 

N 2 5 

\\ 1 9 

S Q R < R .■» > / 2 * R 

S Q R ( R > :? .' 2 * R 



- B /■ 2 * R 



IMR 



G I N R R Y " 

2 * R ::i ; " + " ; S Q R c; - 



- E .' 2 * fl > ! " 



; S Q R < - 
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> R U M 

Q U B D R R T I C E Q U 
PLEASE ENTER 
VALUE D F fl = 1 
VALUE DF B=2 
VALUE a F C = 1 



R T I D N S D L V E R 
THE FaLLDIjJING 



THE 
-1 



RDDTS PRE 



* * D D M E * * 

> R U M 

Q U R D R fl T I C E Q U fl T I D N S D L V E R 
PLERSE ENTER THE FDLLDI.JIN6 
VRLUE QF R=0 
'v' R L U E OF B = 2 
VRLUE DF C=4 

THAT IS R FIRST 
DEGREE EQUATION 

** HONE #* 



• R U r.1 

QUADRATIC E Q U R T I D N S D L '■/ E R 
PLEASE ENTER THE FOLLOWING 
VRLUE □ F R = 2 
VRLUE D F B = 5 
VRLUE D F C = 2 



THE 



RDDTS RRE 



** DDhJE ** 

■RUN 

Q U fl D R fl T I C E Q U R T I D N S D L V E R 
PLEASE ENTER THE FDLLDI.JING 
VRLUE DF ft=4 
VRLUE D F E = 2 
VRLUE D F C = 4 

THE RDDTS RRE IMflGINRRY 
- 4 + 1 5 . 4 9 1 9 3 3 3 S! * I 
-4 - 15.4 9 19333 8 *I 

** DUUE ** 
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Exponential Equation Solver 

Another general routine for solving a particular type of equation is this one 
to solve for an exponent in an exponential equation. 

Given the values A, B, m, and n, this program will solve for x in any 
exponential equation of the form: 

^mx + n ^ B 

For example, the program will solve any of the following problems: 

1. S'' = 40 

2. 5^^ + ' = 7.6 

3. 17''-^ = 8.12 

4. 11I-2X = 247 

If you were to solve an exponential equation by hand, you would probably 
go through the following steps: 

5^ = 40 

log 5^ = log 40 

X log 5 = log 40 

X = log 40/log 5 

X = 1.6021/.6990 = 2.292 
However, in more generalized form, the solution for x is: 

logB 
X = 'og A 



M-(N/M) 



Note that the program to solve this problem is actually divided into two 
sub-programs. The first is a data loader program. It requires that data be en- 
tered in the following order: A, B, m, and n for each equation to be solved. If 
a coefficient is not present, it must be entered as a zero. 

The data for the four problems listed above were entered into the two 
data statements (50 and 60). 

The program as it is presented here can be improved in several ways. First, 
it always solves the same four equations. How can you generalize it to solve 
for other equations? Second, if you were to make use of this routine in an- 
other program, you would probably not be able to use a read statement; how 
could you get rid of it?. 
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10 

DN 
15 

2 

3 
3 5 
40 
M? 
45 

5 

6 
- 2 

99 



PR 

::.D 
PR 
X " 
RE 
X = 
v = 

PR 
TRB 
150 
Dfi 
Dfi 
, 1 
EN 



INT " EXPONENT I RL E Q U fl T I 

LVER " 

I N T " PI B M U 

y :j / M - c N . ' ti > 
^ B <: 1 ) ; 

TO 2 

T R 5 - 4 , 1 r , 5 , 7 . e. , 3 , 1 

TR 17.3.12,1,-3.11,24?', 

D 



>RIJ 
EX 
fl 

5 
5 
1 



N 
PDNENTIRL 

B ri 

4 1 



7 8.12 

1 
1 247 -2 

DRTfl ERROR 



E Q U fl T I D N 5; D L '■/■ E R 
N X 

2.29 

1 .09 



-3 
1 

I N 20 



3.74 
-.65 
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Roots of Any Function 

This program will find the roots of a function, any function! The function 
may be linear, quadratic, cubic, trigonometric or any combination as long at 
it can be represented in the Basic language. The program as it appears here 
finds the roots between -20 and 20 although you can change these boundaries 
in Statement 120. 

The method used involves evaluating the function at small incremental 
intervals, finding places where the value of the function changes sign and 
then, by successive approximations, finding the zero point. This approach 
borrows from Newton's method in the final narrowing down but, unlike 
Newton's method, will not fail to converge in the event one makes an unlucky 
first guess. 

Before running the program you must first type in your function in State- 
ment 100. For example, 

DEF FNA (X) =2*XT3+11*XT2-31*X-180 

DEF FNA (X) = X-4 

DEF FNA (X) =SIN(X)-.5 

You may have to refer to the Basic manual with your system to see exactly 
how a function should be stated. 

The routine used in this program is very powerful and could possibly be 
used as a subroutine in many other programs. 

How can you use this program to help you solve this problem for x? 



I I I r-r= 

x= ^12 +/12 +/12 +/12 yi2 + 
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100 

31 * 

1 10 
120 
130 
+ 1 > 
140 
150 
160 
170 
180 
190 
20 
210 
220 
230 
240 
250 
HEN 
260 
0+. 
270 
2S0 
290 
30 
99 9 



DEF 
X - 1 8 
21 =- 
FOR 
IF S 
MHEN 
K. = I 
J=I + 
IF F 
Z=K 
K = J 
J = Z 
Z = < K 
IF F 
.J = Z 
GDTD 
K = Z 
IF fl 
2 00 
Z = SG 
5 ) .-' 1 
IF Z 
PRI N 
Z1=Z 
NEXT 
END 



F N fi '; K ) = 2 * K ■'^ 3 + 1 1 * K - 2 - 

20 

I=-l 9. 9 TO 20 
G N C F N fl ^: I > > = S S N C F N R •: I 
300 



MR ■:: K ) < FNR < J > THEN 200 



+ .J)y"2 

NR';Z)<0 THEN 240 

250 

B S C F N fl t; Z > ::■ > . 5 T 

N ( Z > * I N T ■: fl B s <: Z > * 1 

000 

=Z1 THEN 30 

T " F >; " ; 2 ; " 5 = " 



■ R U N 

F •: - 5 ■>=o 
F 'C - 4 . 5 :■> 

F-: 4 >=0 
* * B D N E 



•100 
•RUN 
F-: 4 



BEF 



= 

* * 

F H R < X > = X - 4 



= 



** DOME ** 



•100 

•RUN 

F <- 



D E F F N fl <; X ) = S I r-< ''s X ) - . 5 



IS 
1 



-1 



•J 



3 2 
. 23 
.04 
948 
759 
6 65 
2 36 
618 
8 06 
90 1 
. 09 
.18 
. 3 7 



1 6 

2 8 
4 ) 
6 > 
2 > 

> = 

8 :■> 

2 ^ 2 

44~ 
32 



= 

:-' = o 

> = 
= 

= 
= 



= 
= 


■> =0 

> =0 



•* * n D N E * * 
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Plot Any Function 

Here is a nifty program that will produce a plot of any function on the 
screen or printer. 

Before running this program, you must type in your function in line 200. 
Like the previous program, which finds the roots of any function, this one 
will plot any function. You must tell the program between what values you 
want the function plotted, i.e., a minimum and maximum value of the x co- 
ordinate. You also input the x plotting increment you wish. 

As it appears here, the program does not allow the user to select the y co- 
ordinates; the program plots y values between -30 and 30. 

Here are several functions you might want to try plotting: 



Function 

DEF FNA(X) = 2*X 

DEF FNA(X) == 30*SIN(X) 

DEF FNA(X) = X-X*X 

DEFFNA(X) = 30*EXP( -X*X/100) 

DEF FNA(X) = X*X-X 

DEF FNA(X) = Xt2-X-15 



X Limits 


Increment 


-15 15 


1 


-5 5 


.25 


-5 6 


.5 


-30 30 


1 .5 


-5 6 


1 



The last function listed is the one plotted in the sample run with the pro- 
gram. Exponential functions are a great deal of fun and sometimes lead to 
unexpected and interesting results, particularly when combined with trigo- 
nometric functions. Experiment! Have fun! 
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2 
21 
2 2 

23 



25 
29 
30 
3 1 
6 
3 2 



3 3 
3 4 

■ H 

7 

35 
36 



BE 
PR 
I N 
XI 
IN 



IN 
: S 

PR 
PR 
Fa 
IF 



F FNfl 
INT 

F' LI T " 

P LI T " 

P U T " 



INT 
INT 
R X = X 
R B S < 



T 



T = 

IF 
Y = 
IF 
PR 

GO 
PR 



3 7 

3 8 

3 9 

4 1 

42 
43 

99 



PR 
+ 20 



X = 
60 

PR 

NE 

PR 
EN 



FURC. X 
't' < 2 8 

X. I 

Y > 1 4 
INT T 

TD 42 
IHT " 
-- +Y 

INT " 
+ r^o " 
:<-3 
TD 4 2 
INT T 



'; X ) = 
INI T 
F INR 
I NCR 



B-B<. 1 
1 TD 

X > < . 

:j y 2 + 
THE 

THE 
R E C '■«' 



X ^- 2 - X 
IfiL V 
L VRL 
E M E H T 



-15 

n L. U E a 
UE OF 
RL STE 



F 
P 



X =s 



X2 S 
G 

14 

N 3 30 

N 4 10 
^ ; " * " 



= " J X I 
T E P - S 
1 T H E f J 



T fl B ( 1 4 ) 




-Y 

It 

-30 





RB O 



2 - 1 Ci +1 



XT X 
INT T 
D 



4 > 5 " I 
fl B t: 1 2 > 5 '■ X 



" ; T fl E c 
= " ; X 2 



IN 
F I 
IN 



ITIRL VALUE DF X=? 
NRL VlflLUE DF X = -e. 
CREMENTRL STEP =1 

X= 7 



•Y ■ 
■30 



2 - 1 
* 



+10 +20 +30 



+ Y 



X = -6 



** DDHE ** 
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"I'm in the kitchen, dear— using the computer. 
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Sets and Repetitive Trials 



For solving relatively simple problems, the computer may not be any help 
at all. In fact, it may take more time to write a program to solve a problem 
than it would to solve it by hand or with a calculator. This chapter should 
help you recognize problems that are suitable for computer solution and 
those that are not. 

In some problems, you may think that the computer will not be any help. 
However, one thing writing a computer program will always do is force you 
to reason out the approach to solving the problem logically and precisely. The 
computer can't solve problems unless it is told exactly how to proceed; hence 
you must understand a problem completely before you can program it for the 
computer. 

Several of the sections in this chapter discuss sets of data. While the sets 
used in the examples have relatively few elements or values, you should bear 
in mind that real world problems often have thousands or millions of pieces 
of data, and the only practical way to solve problems of this size is with a 
computer. For example, consider how you would most efficiently schedule 
the shipments on a railroad train leaving Boston with 4000 diverse cargoes 
bound for Phoenix and 780 points in between. Now consider that there are 
200 freight trains per day leaving Boston. Now add to that the 10,000 other 
trains leaving other cities every day that must use the same network of track 
and you can see that dealing with real world sets of data is no easy task. 
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Group of Girls and Boys 

In the previous chapter we said that there may be some problems for which 
brute force trial and error is appropriate. This would be the case if the prob- 
lems were relatively small and trying every possible solution would not tie up 
a great deal of valuable computer time. Here is a problem involving two lin- 
ear equations that lends itself to a trial and error approach. 

The problem is as follows: When 15 girls leave a group of boys and girls, 
there are two boys for every girl (lucky girls). Next, 45 boys decide to leave; 
then there are 5 girls for every boy (lucky boys!). How many girls were there 
in the group before anyone left? 

Before rushing to the computer, you must recognize that this problem re- 
quires the solution of two simultaneous equations. If G equals the original 
number of girls and B the original number of boys, then the two equations 
are: 

(G - 15) X 2 = B 

(B-45)x5 = (G- 15) 



10 1=0 

20 FDR G=l TD 100 

30 FDR B=l TD 100 

40 1=1+1 

5 IF 2 * t; G - 1 5 X: > B THEN 90 

6 IF 5 * < B - 4 5 :j <: >< G - 1 5 > T H E N 3 


? PRINT G 5 " G I R L S » " ? B ? " B D Y S 

8 15 D T n 12 

•5 HEY^J B 

10 N E ><; T G 

110 PRHJT "NO SDLUTID^J BET WE 

E tJ 1 8* 10 0" 

12 P R I fJ T I ; " C D M B I N R T I D N S T R 

lED " 

999 END 



■RUN 

4 GIRLS, 50 B D Y ■:. 

39 5 CDMEIfilflTiaNS TRIED 

** DONE ** 



The computer program uses two FOR loops (Statements 20-100 and 30-90) 
to try every combination of values for B and G between 1 and 100 until a 
solution is found or until the program runs out of values. The variable I 
(Statement 40) is a counter which records the number of trials required to 
reach a solution. 
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The program is straightforward and finds a solution after 3950 trials. How- 
ever, it would have been a simple matter to substitute the value of B from the 
first equation in the second one and quickly solve the problem by hand or 
with the aid of a calculator. It is important to recognize that if a problem can 
easily be solved by other methods, (he computer offers little or no advantage. 

Try this problem. You may or may not want to use your computer. If Mat- 
thew can beat Jeff by one-tenth of a mile in a two-mile race and Jeff can beat 
Steven by one-fifth of a mile in a two-mile race, by what distance could Mat- 
thew beat Steven in a two-mile race? (Hint: the answer is not 3/10 mile.) 








#1 X^-^r :!^4^l^f 






''^ 
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Brown's Books 

The use of a trial and error approach can generally be improved signifi- 
cantly if the combinations to be tried can be narrowed down in some way. 
The solution to this problem illustrates how the speed of obtaining a solution 
can be improved well over 100 fold by combining equations and eliminating 
certain solution possibilities. 

Here is the problem. Brown sold 48 books at a fiea market, some for $3 
each, some for $5 each and others for $8. He collected a total of $175. He 
remembered having an even number of S5 books. Can you determine how 
many of each kind of book he had? 

The equations for solution are (letting T equal the number of $3 books, F 
the number of $5 books, and E the number of $8 books): 

T + F + E = 48 

3*T + 5*F + 8*E = 175 

The first program was written simply to try all possible combinations of T, 
F, and E from 1 to 48. It yields three solutions for the problem, although the 
two solutions with an odd number of $5 books can be eliminated leaving just 
the one desired solution. 



10 PRINT "BOnKS PROBLEM" 
2 PRINT "RLL SDL UT I DNS" 

•3 n P D T W T 



2 r,.. .., , 

30 PRINT 

4 PRINT " S 1 " ; T n E ( 1 > ; 
T fl B C 2 > ; " « S " 

5 PRINT 
60 FDR T= 1 TD 43 
7 FDR F=l TO 4 8 
30 FDR E=l TO 48 

? Q „ I F ^ ': T J F + E )< > 4 8 THEN 120 
10 IF t. T * 3 + F * 5 + E * S ) < )■ 1 7 5 T H 
EN 12 
110 PRINT T;TflB(:iO>;F;TflE':20 

12 NEXT E 

13 NEXT F 

14 NEXT T 
9 9 9 END 



BDDKS PRDELEM 

FILL snLUTinr-^s 

« 1 4 5 $ 8 

40 3 S 

37 8 3 

3 4 13 1 



** DOME ** 
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This program took approximately 22 minutes and 1 3 seconds to run on the 
TI 99/2 computer. A typical minicomputer (PDP-8/e) could run this prob- 
lem in about 7.3 seconds. In either case, this is a long time to tie up the 
computer. 

It is rather easy to combine the two equations into one by solving for T. 
The single equation is then: 

2*F +5*E = 31 

In this equation, the Hmits can be reduced (from 48 used in the first run) 
since F cannot possibly be greater than 31/2 or 15.5 and E cannot be greater 
than 31/5 or 6.2. Making the appropriate program modifications leads to the 
second program. 

>6C) 

>70 FDR F=l ID 15 

>30 FOR E=l TD 6 

>90 IF <:F*2 + E*5;i < >31 THEN 120 

MOO 

>105 T=48-F-E 

>140 

> 

>RUN 

BDDKS PRDBLEH 

ALL SOLUTIONS 

* 1 S 5 S 8 

40 3 5 



J7 8 

34 13 

** DONE ** 



? 



Using this program produces a dramatic improvement in the time to solu- 
tion. On the TI 99/2, the time is approximately 2.5 seconds and on the PDP- 
8, about 0.16 seconds. 

Since the problem states that F must be even, a final modification which 
steps F by two in Statement 20, can be made. This version of the program 
takes only 1.25 seconds to run on the Timex and 0.06 seconds to run on the 
PDP-8. 

:>20 PRINT "SOLUTIONS, EVEN Vfi 

LUEDFF" ^^^ 

>70 FDR F=2 TD 14 STEP 2 

Notice the enormous improvement in computing time required for a solu- 
tion, over 1000 fold on the TI 99/2 and 100 fold on a PDP-8. Brute force 
certainly is inefficient! It is generally worthwhile to think through most prob- 
lems, particularly big ones, before rushing to the computer. The computer 
may be fast, but we just improved its performance by 1000 times by using a 
little common sense. 
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Intersection of Sets 

Two sets of numbers can be combined to yield a third set by the operation 
of intersection. The intersection of two sets A and B is the set that contains all 
elements that belong to both A and B. It does not contain any other elements. 
The intersection is usually written A n B. 

For example if M = < 0,2,4,6 > and K = < 1,2,3,4 > , then M n K = 
< 2,4 >. 

This program finds the intersection of two sets of numbers. It has been 
written to find the intersection of the two repetitive sets described in State- 
ments 30 and 40. In the sample run. Statement 30 describes the set 
X = < 1,3,5, ... 19 > and Statement 40 describes the set 
y = < 2,5,8,. . .29 >. 
Notice that successive values of x increase by 2 and y by 3. 

10 PRINT "THE INTERSECTION D 

F SETS" 

2 PRINT " X AND V IS: " 

3 FDR y.= l T a 1 •? STEP 2 

4 FOR Y=2 TO 2 9 STEP 3 

5 IF X = Y THEN 100 

6 NEXT V 

7 NEXT X 

8 :=; T a P 

10 PRINT X 

110 GOTO 70 

9 9 9 E n D 
> 

> R U H 
THE INTERSECTION OF SETS 
X R N II Y IS: 
5 



IJ 



* * D fJ E * * 

However, if the set cannot be so neatly described, it may be desirable to 
rewrite the program to examine any set of data. This is done with the read 
statement which reads into x the data points in the data statement. The pro- 
gram is set up to use the same y set as the first program, but the x set is de- 
fined in the data statement. 

You should be able to see from the first combination of sets that if there is a 
numerical pattern in the sets which intersect, then there is also a pattern in 
the resulting intersecting set. In the example, the x values increase by 2 and 
the y values by 3, hence the values in the intersecting set increase by 2 X 3 = 
6. Although the intersection of these sets could easily be calculated by hand, 
the computer can be an aid in evaluating more complicated sets. 
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10 PRINT "THE INTERSECTION 

F SETS" 

2 PRINT "X R N D Y IS:" 

30 RERII X 

40 FDR Y=2 TD 2 9 STEP 3 

5 ID IF 'f< = V THEN 100 

6 N E >: T V 
70 GDTD 30 
10 PRINT X 
110 i5 D T D 3 

200 DflTfi 2, 3> 8, 9 , 14, 15, 20, ; 
» 2 6 . 2 ? 
99 9 END 



•RU 
TH 
X 

2 
S 
1 



E INTERSECTION 
AND 'T' IS: 



DF SETS 



* DflTR ERRDR IN 30 
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Prime Factors 

A prime factor is a positive integer that has no factor except itself and one. 
The first ten prime factors (or numbers) are 2, 3, 5, 7, 1 1, 13, 17, 19, 23, and 
29. The definition gives the basic method for determining whether a number 
is prime: divide by all smaller integers down to 2, testing whether the remain- 
der is zero for at least one of them. If not, the number is prime. 

But this is highly inefficient. It is obvious that a number is not prime if it is 
any even number greater than 2; hence only odd divisors need to be tried. 
Also, it is not necessary to try divisors greater than the square root of the 
number. 

Since the division method is inefficient, various schemes have been devised 
to avoid division. The basic idea underlying all such schemes is called the 
sieve of Eratosthenes (276 B.C. -195 B.C.). Imagine a list of odd numbers 
from 3 up. Strike out every third number after 3, every fifth number after 5, 
and so on. This will leave only prime numbers. 



10 DIN fl < 1 40 > 

2 C = 

3 PRINT "PRDGRfiri COMPUTES P 
RIME " 

4 PRINT "FRCTDRS DF ANY IHl 
EGER " 

50 PRINT 

6 PRINT "ENTER CZERO;' TD 

5 T D P " 

90 PRINT 

100 INPUT "VDUR NUMBER = " : M 

110 X = 

12 M = I N T < ri :■ 

13 IF M > THEN 160 
13 5 IF M = THEN 330 

140 PRINT "MUST BE POSITIVE" 

15 15 a T D 10 

16 PRINT M ; 
1 ?0 1=1 

180 1=1+1 

190 IF I^'M THEN 250 

2 IF M ^ I< :> I N T < M ^' I :> T H E N ISO 

2 10 X = X + 1 

2 2 fi <: y, ) = I 

2 3 M = M ^ I 

2 4 GOTO 2 ID 

250 IF X=I THEN 300 

2 6 FDR L = 1 TO X 

270 PRINT R<L>; 

2 8 N E >; T L 
290 60TD 3 10 

3 PRINT "IS PRIME" 
310 PRINT 

320 GOTO 90 
3 3 STOP 
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PRDGRRM CDMPIJTES PRIME 
FflCTGRS OF flNV INTEGER 

ENTER ( Z E R a ;^ T D S T Q F' 

V D U R NUMBER =105 
10 5 3 5 7 

'(' D U R NUMBER =72 
72 2 2 2 3 3 

Y D U R NUMBER =89 
S9 IS PRIME 



YDUR NUMBER =4 7 
4 7 IS PRIME 



YDUR N LI M B E R =0 
** DONE ** 



The program here finds the prime factors of any integer, or prints out "N is 
prime" if the integer has no proper divisors. 

Run this program for a large number of different integers and see if you 
can discover relationships between numbers and their prime factors. You 
should also try to figure out the method employed in the program to find the 
prime factors of any integer. To do this, you might want to draw a flowchart 
to show what is happening in the program. This will help you see the method 
used to find a prime factor and might help you in writing a program to gen- 
erate primes. 

In writing a program to generate prime factors, you can use the sieve 
method. However, as the numbers become very large, you will have to figure 
out a way to represent integers with more digits than your computer can han- 
dle at one time. (One approach is described on pp. 19-21 of Computers in 
Mathematics. ) 

Goldbach was a mathematician who made a conjecture that every even 
number greater than 4 can be written as the sum of two prime numbers (16 = 
11 + 5, 30 = 17 + 13, etc.). No one has ever proved it but no one has dis- 
proved ii either. That is why it is called a conjecture. Can you write a pro- 
gram that will prove or disprove this conjecture? Or how about writing a 
program to prove Goldbach's conjecture for even numbers up to 50? You 
should be able to write this program with 12 or fewer statements. 

Here is another problem involving prime numbers. Assume a life span of 
80 years. In what year of the 20th century (1900-1999) would a person have 
to be born to have the maximum number of birthdays occurring in prime 
years? The minimum number? 
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Greatest Common Divisor 

The greatest common divisor of a set of numbers is, as its name implies, the 
greatest integer that will divide into a set of two or more numbers. For exam- 
ple, the set of numbers 12, 20, and 28 have a greatest common divisor of 4. 
Nothing larger than 4 will divide evenly into all three numbers. 

This program will find the greatest common divisor for any set of integers. 
To run it, you simply input the number of integers in your set, type them in 
when requested and let the program calculate the GCD. The heart of the 
calculation is in Statement 1 80. 

Do you know the meaning of a relatively prime set of numbers? Can you 
figure out the meaning from the third sample run of the program or from 
runs of your own? How is a set of relatively prime numbers different from a 
set of prime factors? Can you find a set of 10 integers that is relatively prime? 



10 PRINT "PRDGRftM COMPUTES S 

R E R T E S T " 

20 PRINT "CaMMDN DIVISOR." 

4 INpijT "NUMBERS IN SET ="! 

N 

50 N=INTi::N) 

6 ni N ;< c: 1 > 

70 PRINT "ENTER NUf'IBERS" 
80 :E;=1E25 

9 FDR K = i TO N 

10 INPUT X ■:: K > 

110 IF X c; K ) > S THEN 1 3 U 

12 S = X >: K > 

13 NEXT K. 

14 PRINT 

15 G = 

160 FOR M=2 TO S 

170 FDR 1=1 TO N 

18 IF J< < I .W M < > I N T c; X C I > -' M ) T H 

EN 2 10 

190 NEXT I 

2 6 = M 

2 10 NEXT N 

2 2 P R 1 1^ T " tJ U rf B E R S "5 

230 IF G::>0 THEN 26 

2 40 PRINT "fiRE RELRTIVELY PR 

IM E " 

250 earo 2 70 

260 PRINT "5. CD. =";G 
2 7 P R I t-J T 

2 30 PRirn "ANOTHER SET <Y OR 
N .■^ " ; 

290 INPUT R5 

3 IF fl S = " Y " THEN 30 
310 STOP 

40 



> R LI N 
PRDGRRM CDriPUTES 
COMMON DIVISOR. 



GREATEST 



NUMBERS IN SET 
E r-J T E R N U T1 B E R S 
^ 1 ' 






9 6 



NUMBERS 
fi N D T H E R 



G . C 
SET 



D. 
<■. Y 



^J UMBERS IN SET 
ENTER NUriBERS 
r> 2 
? 36 
? 97 



= 3 



= 1 2 
□ R N > ? 
= 3 



NUMBERS R R E 
RNDTHER SET 



R E L fl T IV E L 'T' PRIME 
f. Y DR N>? ■ 



In the last section we discussed prime numbers. Here is an interesting chal- 
lenge for you involving prime numbers. Until late 1982, the longest pro- 
gression of prime numbers in which all differed by the same number was 17. 
Prof. Paul Pritchard in the computer science department at Cornell Univer- 
sity wrote a program to determine if there was a longer progression. Using a 
DEC VAX-11/780, he found the string of 18 numbers shown below. He also 
discovered fourteen other 17-number progressions and ten 18-number pro- 
gressions, but none yet with 19 numbers. He believes there is at least one; can 
you find it? 



107928278317 


197233324147 


117851061187 


207156107017 


127773844057 


217078889887 


137696626927 


227001672757 


147619409797 


236924455627 


157542192667 


246847238497 


167464975537 


256770021367 


177387758407 


266692804237 


187310541277 


276615587107 
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Cryptarithmetic Problems 

Cryptarithmetic or alphametic problems are arithmetic expressions in 
which the digits are replaced by letters of the alphabet. Each digit is asso- 
ciated with a letter to produce an interesting statement, for example: 

SEND 
+ MORE 

MONEY 

If the college student who sent this message to his father needed $106.52 
for plane fare home, this was the right message to send since this combination 
of letters has one unique solution, in particular: 

9567 
+ 1085 

10652 

However, if the student let things go to the last moment and was in more of 
a rush, he might have reworded the message: 

WIRE 

+ MOR E 

MONEY 

In this case, how much should his dad send? Earlier in the book, trial and 
error approaches to solving problems were discussed. It was noted that the 
brute force approach of trying every alternative was sometimes appropriate. 
Is it in this case? 

No! The number of possible alternative solutions is the factorial of the 
number of different letters in the alphametic expression, i.e., 8! or 40,320. A 
program to try out every one of these possibilities would run for a 1-o-n-g 
time. 

In this case it is much more efficient to apply some common sense to nar- 
row down the number of alternatives. The best approach to this process is to 
divide up the search space into large classes (or sets), according to a common 
property shared by members of each class, and then attempt to eliminate en- 
tire classes by the method of contradiction. 

Consider the "WIRE + MORE = MONEY" problem. Can E = 0? Since 
E + E = Y, Y must also equal 0, contradicting the fact that Y and E must be 
different digits. Thus, the entire class of solutions in which E = can be 
ruled out. 

Consider E = 3. Now Y = 6 and there is no carry to the next column. So 
in this column R + R = EorE+ 10 if a carry is involved. But in either case 
E must be an even number since 2R is always even; this contradicts the 
assumption that E = 3. 

By following this type of classificatory contradiction process for each of the 
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digits in the order E, R, I, O and N, the computer program will search out all 
possible solutions to the problem. Unlike the "send more money" problem, 
the "wire more money" problem has five possible solutions. A smart father 
would choose the lowest solution and wire his son $103.48. 



10 PRINT " P R D G R R t1 FINDS R L L 
S D L U T I a N S " 

2 PRINT "TO THE FDLLDWING F' 
R E L E M : " 

3 PRINT 

4 PRINT " hi 1 R E" 

5 P R I N T " + M D R E " 

6 PRINT " " 

7 P R I N T " M D N E Y " 
SO PRINT 

9 M = 1 

100 FDR E=2 TD 3 

110Y=E+E 

120 IF 'rMO THEN 150 

13 C 1 = 

14 6DTn 170 

15 C 1 = 1 

16 Y = Y - 1 

170 FDR R==l Tp^9^,. 

18 IF R = n THEN 470 

182 IF R=E THEN 47 

184 IF R=Y THEN 47 

190 IF R+R+C1=E THEN 2 20 

200 IF R+R+C1=E+10 THEN 240 

2 10 GDTD 470 

22 C2=0 

2 30 GDTD 2 50 

2 4 C 2 = 1 

25 FOR 1=0 TD 9 

2 6 IF I = M THEN 460 
262 IF I=E THEN 46 
264 IF I=Y THEhJ 46 

26 6 IF I=R THEN 460 

27 FDR D=0 TD 9 

2 3 IF D = t^ T H E f J 4 5 
2 3 2 IF □ = E T H E r-J 4 5 
2 8 4 IF n = Y THEN 450 
2 3 6 IF a = R THEN 4 5 U 
28S IF n=I THEr^ 45 LI 

2 9 N =-- I + a + C 2 

300 IF N>=10 THEN 330 

3 10 C 3 = 

3 2 G □ T D 3 5 

3 3 C 3 = 1 

340 N=N-10 

350 IF N=M THEN 450 

360 FOR W=0 TD 9 

3 6 2 IF W = M T H E rj 4 4 

3 64 IF W=E THEN 44 

3 66 IF W=Y THEN 44 

36 8 IF U=R THEN 440 
:3 7 IF U = r THEN 440 
3 7 2 IF IJ = D THEN 440 

37 4 IF U=N THEN 440 

3 3 IF D + 1 < > IJ + M + C 3 THEN 440 
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3 9 


PR H T '■ 


" ; w ; I ; R ; E 


40 


PR H T " 


" ; M ; n ; R ; E 


4 10 
420 


P P W T " - 




PR NT M ; 


D ; N ; E ! Y 


430 


PR NT 




44 


NEX-- W 




450 


NEX" D 




460 


NEX'-- I 




470 


NEX" R 




413 


NEX" E 
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STOP 





PRDGRflfl FINDS RLL SDLUTIQNS 
TD THE FDLLaUJING PROBLEM: 





W 


I 


R 


E 


+ 


M 


D 


R 


E 


M 


n 


N 


E 


V 




9 


7 


6 


2 




1 


ID 


6 


2 


1 





8 


2 


4 




<a 


2 


7 


4 




1 





7 


4 


1 





3 


4 


S 




9 


5 


7 


4 




1 





7 


4 


1 





6 


4 


8 




9 


2 


8 


7 




1 


V 


s 


7 


1 





3 


7 


4 




9 


5 


8 


7 




1 





8 


7 



10 6' 
* * n D N E * * 



There are other approaches to solving cryptarithmetic problems, but all of 
them benefit greatly from reducing the search space as much as possible be- 
fore putting the problem on the computer. See if you can devise another 
successful approach and write a program to implement it. 

44 



Here are some problems for you to try. 

DONALD TWO 

+ GERALD X TWO 



ABCDE 

X 4 



ROBERT 

ABC 
X DE 

EEC 
DEC 


THREE 

THE 
EARTH 

VENUS 

SATURN 

+ URANUS 


EDCBA 

SPRING 

RAINS 

BRING 

+ GREEN 




PLAINS 



NEPTUNE 

VIOLIN + VIOLIN + VIOLA + CELLO = QUARTET 
THREE + NINE = EIGHT + FOUR 



ONE 

TWO 


FORTY 

TEN 
+ TEN 
SIXTY 


FIVE 
-FOUR 


+ FIVE 
EIGHT 


ONE 
+ ONE 



TWO 
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Sailors and Monkey Problem 

There are many variations of the sailors and monkey problem. Here is one 
of them. 

Five sailors and a monkey were on an island. One evening the sailors 
rounded up all the coconuts they could find and put them in a large pile. Be- 
ing exhausted from working so hard, they decided to wait and divide them up 
equally in the morning. During the night, a sailor awoke and separated the 
nuts into five equal piles, but had one nut left over which he gave to the mon- 
key. He took one pile, hid it, and pushed the other four together and went 
back to sleep. He was followed in this action by the other four sailors, each of 
whom did exactly the same thing. Next morning the remaining nuts were di- 
vided equally with one remaining nut going to the monkey. What is the small- 
est number of coconuts with which they could have begun? 

Although there is an elegant algebraic solution to this problem, a more 
suitable approach for the computer is that of working backwards. A typical 
solution to a problem can be thought of as a path that leads from the given 
information to the goal. However, in this case the goal, or final state, is 
known, thus it is easier to start there and work backwards to the initial state. 

As mentioned at the outset, many sailor and monkey problems exist, in 
fact, an infinite number of them. For example, instead of five sailors, there 
could be three or six or 14. Thus it is desirable to devise a general solution in- 
stead of just one to solve one specific problem. 

In the flowchart and computer program, S is the number of sailors and A is 
the number of coconuts that each sailor received in the final division of the 
pile. Since one coconut was given to the monkey at each division, the total 
number of coconuts left in the morning must be S x A + 1. But this pile came 
from pushing together S - 1 equal piles. Thus, the key condition that must 
hold for (S X A + 1) / (S - 1) to be an integer K, which represents the num- 
ber of coconuts that the last sailor stole from a pile of S x K + 1 coconuts. 
But this pile is the result of pushing together S - 1 equal piles by the previous 
thief, so again (S x K + 1) / (S - 1) is an integer and so on back through all S 
raids on the pile. 

Note in the flowchart (and program) that the first trial value for A is 1 
(Statement 70). In Statement 1 10 this value is increased by 1 until the value of 
(S X K + 1)/ (S - 1) is an integer as tested for in Statement 100. This process 
is then continued until the counter for the second loop, N (nighttime pile di- 
visions) equals the number of sailors. 

Although the program will work for any number of sailors, it takes a fairly 
long time to run for more than five. Remembering what you have learned ear- 
lier in this chapter, can you devise a way to make the program more efficient? 
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10 
LD 
2 
WD 
30 

4 

5 

II ■ 

6 

7 
80 
30 
10 
K + 

1 1 
12 
13 
14 
15 
16 
1? 
TS 
13 
TO 
19 
20 

2 1 
FlC 
2 2 
23 



PP. 

RS 
PR 

RK I 
PR 
PR 
IN 

S 
PR 
R = 
K = 
N = 

I 

1 >-" 
Fl 



G 
N 
K 
I 
G 
P 
TH 

P 
BE 

P 



INT 
AND 
INT 
US" 
INT 
INT 
PUT 

INT 

1 

ft 



F <:. 

CS- 

= n + 

DTD 
= N + 
= CS 
F N 
□ TD 
R IN 
RT" 
RIN 
GIN 
RIN 
RIN 
RIN 
flIL 
RIN 
TDP 



" P R G R R M SDL V E S SRI 
" M a N K. E 'T' P R D E L E ri E V 
"BfiCKtJftRDS , " 
"NUMBER DF SRILnRS= 



S * K + 1 > /• t S - 1 ) = I N T c: .; S * 
1 ;> > THEN 130 
1 
8 

s- 1<: + 1 ) '' <: s - 1 > 

= S THEN 170 

100 
T."THE FEWEST CQCDMU 

t'-SRILDRS can HAVE 

ri 

T " I s " ; s * K + 1 

T 

T "IN THE MDRHING, E 

DR " 

T " G E T S " ■ fi 



>RU 
PR 
MQ 
ER 



N 

DGRR 
N K E Y 
C K W R 



M SQL 
PRDE 
RDS . 



VE3 SfllLDRS RND 
L E M BY IJ D R K I fil 6 



NUMBER DF SRIL.DRS = 3 



TH 
Sfl 
IS 

IN 

GE 



E FE 

ILDR 

7 3 

THE 
TS ? 



WEST 
S CRN 



CDCDNUT3 THAT ; 
HRVE TD BEGIN 



riDRNING, EACH SfllLDR 



* * DONE * * 



>RU 
PR 
MD 
BR 



QGRA 
NKEY 
CKUR 



M SDL 
PRDE 
RDS . 



VES SRILDRS RND 
L E M BY WORKING 



NUMBER DF SAILDRS=; 



TH 
SR 
IS 

IN 
GE 



E FE 

ILDR 

156 

THE 
TS 1 



WEST 
S CRN 
2 1 



c^CDr^lUTS that s 

HRVE TD BEGIN 



MDRNING, EACH SRILDR 



** DDNE ** 
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Super Accuracy 

Under normal circumstances, your computer performs computations to six 
or seven digits of accuracy. Double precision computations increase accuracy 
to 13 digits or so. 

However, it is possible to do computations one digit at a time and assign 
each digit to an element in an array. This will achieve virtually any desired 
accuracy. Any, up to the maximum array size that is. 

This example program performs the rather simple operation of successively 
doubling a number (which is the same as raising 2 to a power). 

If the number to be represented is 8192, then: 

A(4) A(3) A(2) A(l) 

8 19 2 

To add this to itself, first the rightmost digits are added: A(l) + A(l). 
If there is a carry, the variable C is set equal to 1, otherwise C is 0. The total 
is put into B in Line 100. 

If B is less than 10, there is no carry (C=0) and the new A(l) equals B. If 
B is greater than 10 there is a carry (C= 1) and the new A(l) equals B-10. 

This operation is continued for all the digits (D) of the number and, when 
it is finished, the new number is printed in Lines 220-240. 

If A(N) is printed followed by a semicolon(;) for tight packing, the Basic 
print routine would leave a space in front of each digit (for the sign) and a 
space after each digit (for readability). In the program here, these spaces are 
not wanted, hence the print routine in Line 230 is used which prints the string 
value of the ASCII value of each digit (which is the same as the digit itselO 
but without the spaces. 

This program, incidentally solves the challenge to calculate the number of 
moves in the Towers of Brahma problem (see "Change for Any Amount to 
$5.00). The approach is also used in the next section, "Palindromes." 

10 PRINT "COMPUTES 2 TD Nth 

TD ANY" : "flCCURfiCY. " : : 

2 DIM fl < 1 > 

:?. M = 

40 C = 

50 D = 

€.0 R •: 1 ) = 1 

70 1=0 

SO C = 

90 M = ri+i 

100 1=1+1 

1 ID B = i=i<:i > + fli; i>+c 

120 IF B>9 THEN 150 
130 C=0 



140 GQTD 170 
150 B=E-1 



50 



>RUN 








CgMPUTES 2 TO Hth 


TO ANY 




RCCURfllCY. 

H 

A lt= 160 






C = l 




5 32 170 


fi < I > = B 




6 64 180 


IF KB THEN 100 




7 128 lyQ 


IF C=l THEN 100 




S 25 6 200 


PRINT MS 




9 512 210 


D=I 







10 2 4 2 2 


FDR N=I TD 1 STEP - 


1 


1 


2 4 8 2 3 


PRINT C H R S C Fi < M ) + 4 8 > 


5 


2 


4096 240 


NEXT N 




.3 


8192 250 


PRINT 




4 


16384 260 


GOTO 70 






32768 






] g 








7 


13 10 7 2 






18 


262144 






19 


524288 






20 


104S576 








2097152 






22 








23 


3 3 88 6 08 






2 4 


16777216 






25 


3 3554432 






26 


67108864 






2? 


134217728 






28 


268435456 








536870912 






3 


1073741824 






3 1 










4294967296 






33 


1353993459 2 






3 4 


171 79869184 






35 


34359738368 






36 


687194767 3 6 






37 


137438953472 


N,^-"-.''*' ^V^ 




-T2" 


^^15T5 99627370496 
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90 071992 5 4740992 
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18 14398509481984 




55 


36028797018963968 




56 


720575 9 403792 


7936 




57 


14411 518S075855S72 




58 


28823037615171 1744 




59 


57 6 46 75230 3 42 J 48 8 




60 


1 152921504606846976 






2 3058430092136 9 3952 




62 


461 1686018427 


387904 




63 


9 22337203 68547758 OS 




64 


13446744073709551616 




65 


36893483 1474 19103232 




66 


73786 9 76294 8 3 8206464 




67 


1475739525896 


7641 2928 




68 


295147905 1 793 


52825S56 




69 


590295810358705651712 




70 


1 180591620717 


41 1303424 




71 


236HS3 24143482260684S 




72 


47223664?: 28696452136 9 6 
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Palindromes 

A palindrome is a word, verse, or number that reads the same backwards 
or forwards. For example, the words "mom" and "eye" are palindromes. So 
are each of the lines in this verse: 

Egad, a base life defiles a bad age 

Doom an evil deed, liven a mood 

Harass sensuousness, Sarah 

Golf; No, sir, prefer prison-flog 

Ban campus motto, "Bottoms up, MacNab" 

Numeric palindromes are those numbers which read the same backward as 
forward. The examination of these numbers is a field rich with possibilities 
for creative computing. 

One conjecture concerning palindromes raises an interesting unanswered 
question. Begin with any positive integer. If it is not a palindrome, reverse its 
digits and add the two numbers. If the sum is not a palindrome, treat it as the 
original number and continue. The process stops when a palindrome is ob- 
tained. For example, beginning with 78: 

78 
+ 87 



165 
+ 561 

726 
+ 627 

1353 
+ 3531 

4884 

The conjecture, often assumed true, is that this process will always lead to 
a palindrome. And indeed that is just what usually happens. Most numbers 
less than 10,000 will produce a palindrome in less than 24 additions. But 
there is a real thorn in the side of this conjecture, the number 196. Can you 
determine if a palindrome will ever be produced with a starting number of 
196? 

The number 196 will produce 1675 after two reversals, but after 100 rever- 
sals the resultant sum has 47 digits and is still not palindromic. Why mention 
1675? Because ten other numbers under 1000 will also lead to the sum of 
1675 and thus may not become palindromic. The first five of these numbers 
are 196, 295, 394, 493, and 592. What are the other five? 
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The program here will accept any number as a starting value and complete 
the process of adding the successive reversals and testing if the sum is a palin- 
drome. Try it with some numbers and see if you can identify any patterns. 



10 PRINT "PRQGRflM TESTS CDHJ 
E C T IJ R E " 

20 F'Rir^T "THRT REVERSING DIS 
ITS fl N D " 

3 PRIMT "flDDIHG WILL fi L l/J fl Y S 
PRDHUCE" 

4 PRINT "ft PRLINDRDME." 

45 PRINT "ENTER tZERQ) TO 
QUIT " 

5 HIM B ■:: 5 > 

6 PRINT 

7 INPUT " S T R R T NUMBER = " : R 
75 IF fl=0 THEN 440 

7 7 PRINT 

SO E = 

90 E = E+-1 

10 fl = R .■•' 1 

110 IF INT(R>>0 THEN 9 

120 FOR C=E TO 1 STEP -1 

130 fl=R*10 

14 B C C ) = I N T <; fi - 1 * I tJ T ■:; fl .' 1 > > 

150 NEXT C 

16 D = 

170 FOR C=l TD INT<E-'2> 

18 IF E '' C > = E < E + 1 - C !) T H E f< 2 

ISO D=l 

2 NEXT C 

210 FDR C=E TD 1 STEP -1 

2 2 P R I r^l T C H R S < E < C ) + 4 8 > 5 

2 3 H E X T C 

240 IF D=l THEhJ 270 

250 PRINT " IS R PRLINDRQME" 

26 GDTD 6 

270 PRINT " IS NOT R PflLINDR 

D ri E " 

2 3 IF E /' 2 < = I N T < E -^ 2 ) T H E N 3 
2 9 B C 1 1 J T t E -' 2 > + 1 :> = 2 * B C I N T >; E .' 

2 > + 1 > 

30 FDR C=l TD INT(E.'2) 

3 10 B C C :? = B C C > + B C E + 1 - C ) 
3 2 NEXT C 

3 3 FDR C = 1 T D ir^l T C E ,' 2 ) 
3 4 B C E + 1 - C > = B O: > 

3 5 NEXT C: 

o' 6 E C E + 1 > = 

37 FQR C=l TD E 

8 E <: C + 1 :> = B <: C + 1 > + I N T c: B < C: ) -^ I 

9 B <: C > = B < C > - 1 * I N T <: B ■; C ) .'■ 1 

4 00 NEXT C 

4-^0 E = E+1^"^^ ''^"""^ THEN 160 
4 30 GDTD 160 
4 4 END 



fl 



) 
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> R IJ N 
PRDGRfiM TESTS CONJECTURE 
THAT REVERSING DIGITS flUD 
flDDINe WILL RLUfl'T'S PRDBUCE 
fl PRLIhJDRDME. 
ENTER •;: Z E R D ) T D QUIT 

START NUMBER =19 

19 IS UDl Pi PALI N drome 
110 IS NOT fl P fl L I r J D R D M E 
121 IS fl PALINDROME 

START NUMBER =96 

9 6 IS NOT A PALINDROME 
165 IS NDT fl PflLINDRDME 
7 2 6 IS NDT A F' A L I tJ D R M E 
1353 IS NOT A PRLINDRQtIE 
4 8 84 IS A PALIf^tDFlOME 

START NUMBER =0 

Using this method, write a program that examines all the integers between 
1 and 10,000 excluding those that sum to 1675 at any point. What does this 
show? By the way, you will have to devise a way to deal with 14-digit integers 
which are larger than your computer can normally handle. 

Huh? Is this program "too hot to hoot?" 
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Convergence and Recursion 



The computer is especially suitable for doing repetitive and tedious calcula- 
tions. Two mathematical approaches for solving problems that involve repet- 
itive calculations are convergence and recursion. 

Some problems can be reasonably easily stated in words or described with a 
few simple equations but there are many possible solutions. For example, how 
many ways can you make change for a dime? It is simply stated and the num- 
ber of ways can be enumerated fairly easily: two nickels, one nickel and five 
pennies, or ten pennies, three ways in all. But if you want to solve for all the 
ways of making change for a dollar or five dollars, it would be nice to have 
some help. 

Help on this kind of problem comes from a class of computer program that 
simply breaks the problem into smaller ones and counts up all the alternative 
solutions according to a set of rules. But an even more powerful technique is 
known as recursion. Using this technique, a simple solving algorithm or rou- 
tine is set up to solve the smallest subset of the problem. The unique power in 
a recursive routine comes from the ability of the routine actuallly being able 
to call itself. This is discussed further in the second program in this section. 

Another approach for solving problems that do not have an exact answer is 
that of successive approximations. For example, the exact value of pi, e or the 
length of an irregular curve cannot be precisely determined. But by means of 
increasingly accurate approximations, it is possible to approach the desired 
value from above or below or to converge on it from two directions. The last 
four programs in this chapter illustrate successive approximations and 
convergence. 

55 



Change For a Dollar 

Even though there is very little you can buy for a penny these days, the 
coin will probably be around for some time to come since it is needed to make 
change for odd amounts of sales tax and to fill up penny collections. 

Today U.S. coinage consists of five coins: penny, nickel, dime, quarter, and 
half dollar. How many ways can coins of these denominations be used to 
make change for one dollar? For example, one way is two half dollars, an- 
other is one half dollar and two quarters, and so on. Make a best guess now 
and write it down before you read further. 

There are several different ways to approach a problem of this kind. One is 
to break it down into smaller, more easily solved problems. In other words, 
how many ways can you make change for a quarter? For a dime? You would 
solve these subproblems and combine the answers to give the overall solution. 

If you were more mathematically inclined, you could write a series of equa- 
tions relating each piece of change to every other one and to the dollar and 
solve them. 

A third approach is to do the problem by writing down combinations until 
all the different possibilities are exhausted (or until you are exhausted) and 
then count them all up. This might be called solving the problem by exhaus- 
tion and is a method quite suitable for putting on the computer. 

Write a program that uses this approach to solve the problem. If you use 
loops and count by one, it could take a long time for the computer to run 
through all the possible combinations, possibly many hours. 

Also, if you want to print out all the possible combinations, be warned that 
the printing could also take quite some time and a fair amount of paper. 
There are more combinations than you might think! 

In fact, most people will not be able to guess the answer to this problem, or 
even come close. Ask several of your friends how many ways they think a 
dollar can be changed. Record all the responses and then tabulate them on 
your computer. What is the mean (average) of all the guesses? The extremes? 

The program included here uses the first method discussed to solve the 
problem, in particular, breaking down the problem into subproblems and 
then combining the solutions into one final answer. 

First, the main problem is broken into the next smaller one of making 
change for half dollars. There are three such problems: no half dollars (H = 
0), one half dollar (H = 1), and two half dollars (H = 2). The last problem is 
trivial since there is only one way, but the other two need to be broken down 
further. 

This is done by dividing the remaining money into quarters and consid- 
ering the subproblems on down to the lower denominations. As the number 
of subproblems is expanded, each one becomes easier to solve. In fact, 
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subgoals, which can be solved in only one way, are finally reached. For exam- 
ple, if H = 1, Q = 1, D = 2, and N = 0, then the pennies (P) must equal 
five in order that the total equal 100. 

Notice that at the quarter, dime and nickel stages, adjustments are made in 
the limits of the loops depending upon how much money there is left to 
change. For example, if H = 1 , the only possible subgoals for quarters are 0, 
1, and 2, but not 3 or 4. Also notice that there is no need to test combinations 
of coins to see if they add up to 100, nor is it necessary to include the penny 
as a variable. Simply counting the number of subgoals is sufficient since each 
one can be solved in only one way. 



10 C = 

20 FOR H=0 TD 2 

3 FDR 0=0 TD 4-2*H 

40 FDR D=0 TO 10-5*H-.5*Q 

50 FDR N = TD 2 - 1 U* H - 5 * Q- 2* 

II 

6 C = C + 1 

7 N E >; T H 

8 NEXT D 

9 NEXT Q 

10 NEXT H 

110 PRINT "THERE flRE"JC!i"DIF 
FERENT WflY£":"TD CHANGE ft DO 

L L Fl R B I L L . " 



>RUN 

THERE ARE 292 DIFFERENT WHYS 
TD CHANGE fi HDLLAR BILL. 

* * n D N E * « 



Try to make some changes in this program or write a new one to solve the 
following problems. Say you want two quarters in your change to play some 
video games. In how many ways can a dollar be changed to provide at least 
two quarters? 

Visiting a small town, you find the parking meters still take pennies. In 
how many ways might you get change so that you had at least three pennies? 
Is this any different than the number of ways that would give you five pen- 
nies? Say you want to make a phone call also; in how many ways can you 
change a dollar to produce at least four pennies and one dime? 
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Change For Any Amount to $5.00 

Another way to attack the change problem in the previous section is by 
means of the programming technique called recursion. Get familiar with this 
one — it is very powerful! Donald Piele and Larry Wood described this 
method in an issue of Creative Computing. 

First, define the variables which represent the number of ways to make 
change for n cents using the coins specified: 
A Only pennies 
B Nickels and pennies 
C Dimes, nickels, and pennies 
D Quarters, dimes, nickels, and pennies 
E Halves, quarters, dimes, nickels, and pennies. 
Initially, there are two subproblems in making change for n cents. In the 
first, no half dollars are used, and D is the number of ways to change n cents. 
Second, when one or more halves are used, after one is paid, there remain 
n-50 cents to pay which can be done in En-50 ways. 

Since these two cases are mutually exclusive, it can be inferred that £„ = 
Dp + En_5o. Similarly, 

Dn = Cn + Dn 25 

Cn = Bn + Cn-10 
Bn = An + Bn_5 

Now, begin with the simplest case and build up to E|oo- First of all, it is 
easy to understand why E,, = 1. From above, when n = 50, E^^ = D^g + 
Eg, and it is possible to make change for 50 cents only one more way if half 
dollars are allowed. Therefore E^ = 1 . Likewise, D^ = C^ = B^ = A^ = 1 . 
It is also true that An = 1 for all values of n since there is only one way to 
make change using only pennies. Now the recursive relationships can be used 
to solve the original problem. 

This is the strategy used in the program. It also has the added advantage 
that it can count the number of ways of making change (with coins) for any 
specified amount. 

Now it is your turn. Can you modify the program here to include dollar 
bills so it could count the number of ways to make change for any amount up 
to $10.00? 

Using any method of change making you prefer, write a program to make 
change for one ruble. Russian coins come in denominations of 1, 2, 3, 5, 10, 
15, 20, 50, and 100 kopecks. There are 100 kopecks in one ruble. 
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10 
flVS 
BM 
20 
3 
40 

5 

6 

7 
SO 

" ! K 

■9 

100 

1 1 

120 

130 

140 

150 

160 

170 

li=.0 

I 90 

200 

210 

220 

230 

240 

250 

260 

270 

280 

29 

6E 

FER 

3 



PRINT 
TO " : 
DUNT" 
PRINT 
HIM fl 
II I M 
D I M 
D I M 
D IM 



PRDGRflM COMPUTES h) 
MAKE CHRNijE FDR RNY 
"' TO S 5 . " 



UP 



B< 
C < 

D < 

E 



INPUT 



M = I N T 



R 

B 

E 



101 > 
101 > 
101 > 
10 1) 
1 1 ) 
" AM DUNT 



J + 1 



TO CHfiNGE = 



FDR 



■J ) 



fl 

B<: 
£<: 

IF 
C < 

11 <■: 

IF 

Ii< 

EC 

I F .J 

E <: J > 

NEXT 

PRINT 

PRINT 

FOR t 

E N T Ul fl Y 

STOP 



J 
J 

J 
J> 

J .5 

.1 :> 
J> 



< 2 * X 
= 1 
= 1 
= 1 
= 1 
= 1 
J = 2 
= 1 
) =fl C J 

) = B ■: J 

" THEN 



Tn M 



= E C 
= C C 
< = 5 

= cc: 

= B< 

< = 1 

= D< 
J 



J-l ) 
210 



J ji + c c J - ; 

^ f H E N 2 4 
J > + D < J ~ 5 > 



THEN 270 
+ E < J - 1 ) 



V a U 



CAN 
I N " 



MRKE CHRN 
E < t1 ) ; " D I F 



>R 
P 
H 
U 



UN 
RQG 
ARE 
P T 



RAM COMPUTES WAYS TD 

CHANGE FDR ANY AMOUNT 
D S 5 . 



AMOUNT TO CHANGE =4 

YDU CAN MAKE CHANGE 
S 4.99 IN 57420 
D I F F E R E tJ T U fl Y S 

* * DONE * * 



99 
FOR 



>R 
P 

M 



UN 

R OG 
AKE 
P T 



RAM COMPUTES JJAYS TO 

CHANGE FDR ANY AMOUNT 
D :E 5 . 



AMaUNT TO CHANGE =1.13 



□ U 
1 



CAN riAKE CHANGE FDR 

13 IN 3 8 4 DIFFERENT WRY: 



* * D D r J E * * 
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Perhaps the most famous problem used to demonstrate the principles of 
recursion is the Towers of Brahma. It is sometimes called the Towers of Ha- 
noi or Pharoah's Needles. Here is the problem in as close to original form as 
possible. You should be able to solve it with a relatively short program using 
recursion. 

In the great temple at Benares beneath the dome which marks the center of 
the world rests a brass plate in which are fixed three diamond needles, each a 
cubit high and as thick as the body of a bee. On one of these needles, at the 
Creation, God placed 64 discs of pure gold, the largest disc resting on the 
brass plate and the others getting smaller and smaller up to the top one. This 
is the Tower of Brahma. 

Day and night unceasingly, the priests transfer the discs from one needle to 
another, according to the fixed and immutable laws of Brahma. These laws 
require that the priest on duty must not move more than one disc at a time 
and that he must place this disc on a needle so there is no smaller disc below 
it. When the 64 discs shall have been thus transferred from the needle which, 
at the Creation, God placed them, to one of the other needles; tower, temple, 
and Brahmans alike will crumble into dust, and with a thunderclap, the 
world will vanish. 

If the priests were to effect one transfer every second, and work 24 hours 
per day for each day of the year, it would take them 58,454,204,609 decades 
plus slightly more than six years to perform the feat, assuming they never 
made a mistake — for one small slip would undo all the work. 

How many transfers are required to fulfill the prophecy? Try out your pro- 
gram with fewer discs than 64 to make sure you are on the right track. Here 
is a table of the first few transfers: 



iscs 


Moves 


Discs 


Moves 


1 


1 


6 


63 


2 


3 


7 


127 


3 


7 


8 


255 


4 


15 


9 


511 


5 


31 


10 


1023 
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Converge on e and Pi 

An incredibly important mathematical constant is designated by the small 
letter e. This constant is both irrational and transcendental. Look up those 
terms in a dictionary or math book if you wish, or just plunge on to the next 
paragraph. 

The constant e was first derived by John Napier, also the inventor of loga- 
rithms, to whom we owe an eternal debt of gratitude. Why? If e had never 
been discovered, advances in mathematics, physics, and astronomy would 
have lagged a century or more, because e is the base of all natural logarithms 
and these logarithms are the basis for many branches of science and 
mathematics. 

How is e calculated? The constant e is the limiting value of this expression 
as n approaches infinity: 

e = (l + Vh)" 

Its exact value can never be found, but to 15 places e equals 
2.718281828459045 . . . How can e be calculated? First take l-'/a and square 
it; that equals 2-'/,. Then cube l-'/j and you get 2.3686. Raising 1-^, to the 
fourth power gives 2.414, and so on. Write a program for this method and 
have it print out the initial value of e and each value after each additional 
fraction is added on. 

Another approach is to expand the expression above using the binominal 
theorem and, agam, letting n approach infinity. The expression for the expan- 
sion is: 

e= 1 + \ + 2! + 3! ■•" 4i + ?! ■ • • ;^ 

Now here is where the computer can again be of some assistance. Since 3! is 
3*2! and 4! is 4*3!, all the calculations need not be done for each additional 
fraction. Look at the program and particularly note the calculations in State- 
ments 70 and 90. 

D N E " : " H I 

2 . 6 6 6 e 6 6 1=. 6 7 
2 . 7 8 3 3 3 :?. 3 3 
2 . 716 €.6 6667 
2. 718055556 
2.71 8253968 
2. 7 1827877 
jv,^, _!..„,.. 2.718281526 

110bDTn40 2.7182818 01 

2 . 7 1 8 2 8 1 S 2 6 
2.718231328 
2 . 7 18 2 8 18 2 8 
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10 PRINT " 


CONVERGE 


T ' F C T N -' 4 ' 


TO 


HRLT 


20 E=l 






30 1=0 






40 1=1+1 












isd FOR J=l 


TD 


I 


7 D=D*J 






80 NEXT J 






90 £ = £+1-^1; 






100 PRINT 


E 





Pi is another important mathematical constant that is irrational (meaning 
its exact value can never be determined) and transcendental (meaning it is not 
the solution to any algebraic equation). Interestingly, pi was known and used 
by the ancients. Archimedes, who lived in the second century B.C., by using a 
regular polygon of 96 sides (nearly a circle), proved that the value of pi was 
less than 22/7 and greater than 3-^%!, a remarkable achievement for the 
mathematics of his day. 

Ptolemy in 150 A.D. used the value of 3.1416 for pi and in the middle of 
the sixteenth century the amazing fraction 355/113 was discovered, giving 
the value of pi accurately to six decimal places. 

Incidentally, in 1 897, the General Assembly of Indiana passed a bill ruling 
that the value of pi was four. 

Several infinite series can be used to grind out increasingly accurate values 
for pi. One such series is (l-Vs + Vs - V? + Vg . . .). This series is called an 
arithmetic series and converges very slowly. The program here displays only 
every 500th value of the series. Don't be alarmed if the program does not 
seem to be running very fast; a fair amount of calculating is going on between 
each value that is printed. 



10 
H R 
20 
S S 

3 
LT" 

4 

5 
60 
70 
SO 
SO 

1 
1 10 
120 
130 
140 
150 



PR I 
N " : 
PR I 
-L- 

PRI 

■ II 

S=l 

1 = 1 

Q=0 

P = 

Q = Q 

P = P 

1 = 

S = 

IF 

Q = 

PR 

GD 



NT "CONVERGE D N PI WIT 

"INFINITE SERIES" 

NT "BE PATIENT; THIS I 

n - bl" 

NT "HIT • F C T N .^ 4 •' T D H fl 



CONVERGE □ N PI WITH R N 
INFINITE SERIES 

BE PRTIENT; THIS IS S-L-D-W 
HIT ' F C T N .-■• 4 • T D HALT 



+ 1 
+ S/ 

1+2 

Q < 4 9 9 THEN 80 


I r^ T P * 4 
TD SO 



.3 


14 3 5 9 6 6 6 


3 


. 14059065 


'J- 


14 2 2 6 6 5 6 


3 


14 10 9 16 5 2 


3 


141993455 


2 


14 12 5 8 6 5 2 


3 


14 18 7 8 9 4 


3 


1 4 1 ::; 4 2 1 5 3 


3 


1 4 1 S 1 5 3 2 1 


3 


1 4 1 :5 9 2 2 5 3 




141774836 


3 


14 14 2 5 6 5 3 


3 


141746808 


r^ 


14 144951 


3 


14 17 2 6 2 5 4 



Actually, the approach used by Archimedes converges much more quickly 
and, with the aid of a computer, it is possible to go far beyond the 96-sided 
polygon used by Archimedes. 
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His approach was to construct inscribed and circumscribed polygons and 
measure the perimeters to approximate the circumference of a circle. Con- 
sider a polygon circumscribed around a circle of radius 1. The perimeter 
equals the length of one side times the number of sides. Since the tangent of x 
= AB/BC, but BC = 1, then tan (x) = AB and the length of a side = 2 tan 
(x). Since the circumference = 2 pi r and r = 1, then pi is the circumference 
(or perimeter of an n-sided polygon) divided by 2. 

Similar trigonometry leads to the perimeter of an inscribed polygon being 
equal to the number of sides times sin (x) * cos (x). 

The second program produces values for pi using inscribed and circum- 
scribed polygons. Unfortunately, there is one large flaw in the program, be- 
cause degrees must be converted into radians in Statement 50. This means, of 
course, that you must already know the value of pi, since the conversion fac- 
tor is 360 degrees divided by 2 pi. 

Setting this flaw aside, it is interesting to note how quickly this program 
converges on the value of pi compared to the preceding one. That is because 
this one converges geometrically rather than arithmetically. 

Can you figure out a geometric convergence to the value of pi that does not 
require that you know it (or a conversion factor) before you start? 

10 PRINT "CONVERGE DN PI H 

IT ' F C T N -' 4 ■' T D H R L T ' 

2 PRINT " I r^ S C R I E E D C I R C U 
M S C R I B E D " 

3 N = 6 

4 N = 2 * N 

50 X=36 0/<N*57.29578> 

6 PRINT N * S irj < X :• » C D S < X ) .•-■ 2 , N 
* T fl N C X > '' 2 

7 G □ T D 4 

CONVERGE D N PI 

HIT '■ F C T N .^ 4 ■' TO HALT 

INSCRIBED CIRCUMSCRIBED 

2.598076 198 3.464 10 15S 

2.9 99999 9 77 3.21539 28 1 

3. 105 S 28515 3.15965 9 915 

3.132628537 3.1460S6188 



139 3 5 0176 3.142714573 
14 10 3 19 2 4 3.141873023 
141452446 3.1416 6 272 
1415 57581 3.141610 15 
3. 1 41597 00 8 
:3 . 14 15 9 3 7 2 2 
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■::! . 1 4 1 .I' 5 1 ' ? » 1 
3 . 14 15 8 3 8 6 5 
3. 1415 9 0437 



Convergence on Pi Revisited 

In answer to the question posed in the last paragraph of the preceding sec- 
tion, here is a way to converge on pi without knowing its value beforehand. 

As in the previous program, the basic approach is to add up the length of 
the sides on an inscribed polygon and divide by 2r to obtain a value for pi. 
The program starts with a square (four sides) and doubles the number of sides 
each time. If the old side length is S, then the length of S' of a side of a new 
polygon with twice as many sides is obtained by applying the Pythagorean 
theorem. In particular, 

X' + (S/2)' = R' 



(R-X)' + (S/2)' = (S')' 



Thus, 



S' = \/(R- VR' - (S/lff + (S/2f 

It is easy to reduce this formula algebraically, but accuracy suffers if this is 
done. Also, you will find that S*S is slightly more accurate than St 2. 

Unfortunately, striving for maximum accuracy is somewhat moot on a 
computer that does not have double precision arithmetic. Notice that ac- 
curacy does not improve with more than8192sides and, indeed, as numbers in 
the calculations start to exceed the capacity of the computer, the accuracy 
starts to deteriorate badly. 

There is yet another method to compute pi by convergence. It uses discov- 
eries of Gregory and Euler. Gregory discovered the formula for arctangent: 
Arctan x = x - xV3 + xV5 - xV7 + ... 

Euler came up with a rather interesting formula for pi: 
pi = 4 (arctan CA) + arctan {%)) 

See if you can combine these two formulas to calculate pi. If you are very 
clever, you can do your calculation to yield far more than the seven decimal 
place accuracy obtained by the programs presented here so far. 
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PRINT "CONVERGE DN PI WIT 
"INSCRIEEIl PDLVGDNS":" " 
PRINT "WHEN PROSRfln BRERK 
T Y P E " : " ' C D N •' TO C D N T I N U E " 



PRI 

R = l 

Z = 2 

N = Z 

S = R 

FDR 

I F 

BR 

PR 

V = 

X = 

s = 

N = 
NE 



NT 



+ z 
* sn 

K = 
K. < > 
EflK 
I NT 
S*S 
R-S 
SQR 
N*2 
XT 



SIDES 



PERIMETER 



R'^Z ) 
1 TO 2 S 
14 THEN 
1 10 



1 10 

N ^C S * N > ' ' '.■ Z * R ) 
^ < Z + Z > 
Q R < R * R - V y 
(. X * X + V -J 



K 



> R LI H 
CONVERGE DN PI WITH 
INSCRIBED PDLYGDNS 



WHEN P R D G R fl \'\ 


B R E Fl H:: S .- TYPE 


■CON* TD CONTINUE 


SIDES 


PERIMETER 


4 


2.8284 2 7125 


s 


3. 0614 6 7459 


41 


3 . 12 14 4 5 15 2 


3. 13 6548491 


64 


3 . 14 3 3 115 7 


123 


3 . 14 12 7 7 2 5 1 


256 


3 . 14 15 13 8 1 


512 


3 . 14 15 7 2 9 4 


1024 


3 . 14 15 8 7 7 2 5 


20 4 8 


3 . 14 15 9 14 2 2 


4096 


3 , 14 15 9 2 3 4 6 


8 192 


3 . 14 15 9 2 5 7 7 


16334 


3 . 14 15 9 2 6 3 4 


* break; POINT 


RT 1 00 


CON 




32768 


3 . 14 15 9 2 6 4 9 


6 5 5 3 6 


3. 1415 9 2652 


13 10 7 2 


3 . 14 15 9 2 6 5 3 


2 62 144 


3 . 14 15 9 2 6 5 4 


52423S 


3. . 14 15 9 2 6 5 4 


104S576 


3. 141592654 


2 9 7 15 2 


3 , 14 15 9 2 6 5 4 


4194304 


3 . 14 15 9 2 6 5 4 


S381360S 


3 . 14 15 9 2 6 5 4 


167772 16 


3 . 14 15 9 2 6 5 4 


3 3 5 5 4 4 3 2 


3. 14 15 9 2 6 5 4 


67108864 


3 . 14 15 9 2 6 5 4 


1 34217728 


3 . 14 15 9 2 6 5 4 


26843545 6 


3.141 592654 


536870912 


3. 141592654 



** DONE ** 
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Length of Any Curve 

The previous programs have demonstrated how it is possible to compute a 
very accurate value of pi by adding together the length of the sides of a poly- 
gon as it approaches a circle and dividing by 2r. 

Using a similar approach, it should be possible to inscribe a polygon, or 
portions of a polygon, inside any regular curve and thus determine the length 
of the curve. This program approximates the length of any curve as defined in 
Statement 110 by dividing it into an increasing number of subintervals and 
computing the sum of the secants (a straight line that cuts a curve at two or 
more points). 

To run the program, you must enter the formula or equation describing 
your curve in Statement 1 10 in the form; 

DEF FNA(X) = Any function of X 

You then type RUN and enter the end points of the curve you want to use 
in your calculation. These points are entered in the form of the abscissa, 
which means the horizontal (or y) coordinate of the point. 

The program is written to sum the successive secant lengths and to cal- 
culate the percent of change in each summation compared to the preceding 
one. The sample run uses the function 2x^ + 3x^ -2x +3. Note the substan- 
tial improvements in the length calculation as the number of intervals in- 
creases from 2 to 16, but the rather slight improvements beyond 16. 

Try this program with different curves and functions. It might help to plot 
the function first (remember the program to do that?) and then compute its 
length. Is this method of length calculation more accurate for a function with 
no changes in direction of the curve within the interval selected or for a func- 
tion with one or more changes? Why? 
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RINT "CURVE LENGTHS"!" " 
PINT "DEFINE YDUR FUNCTI 
N " : " L I N E 110":: , „ ^ ^ . „ 
RINT "ENTER flBSCISSflS OF 
PTS":"<SMRLLER DNE FIRS 



NPUT P , Q 

RINT 

PERCE 

RINT 

GE IN 

RINT 

TH" 

RINT 

1=0 

FOR M 

DEF F 



DF SUM 
S E C R N T 
TERVflLS LENGTHS 



NT" 

" S U E I N - 



□ F 



= 1 TD 9 

N fl <; X > = 2 * X -^ 3 + 3 * X •"• 2 - -^ 



E = 2 '^ C 
H=<Q- 
S = , 
FDR I 

L = S. Q R 
H> > '^2 
S = S + L 
NEXT 
IF SI 
PRINT 
PREV 
GDTD 
P 5 = < < 
PRINT 

S1=S 

N E X T 



N-1 ) 
P>^E 

=0 TD E-i 

f CFNfi<P+I*H + H)-FNR C 

+ H*H) 



>0 THEN 220 

E;TRBt:5:';s;TfiBc;i8> 

V fl L " 
240 
RBS<Sl-S;o.'Sl>*100 

E;TiRB<5;';s;TflB«;i7) 



N 



CURVE LENGTHS 

lUR FUNCTION IN 



DEFIN 
LINE 

ENTER 
CSMRL 

ND . D 
S IJ B I N 
T E R V R 

1 
2 
4 

8 

16 
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E YOl 
1 10 

RBS 
LER 

F S 

S 

L3 L 

525 . 
525. 
529. 
531 . 

53 1 . 

5 32. 

5 32. 



CISSftS 
□ NE FI 

UM aF 
E C fl N T 
ENGTHS 

0466641 
158262 
65224? 
17 13 3 

9 6 4 2 5 5 

16 5 6 5 

041678 



DF END 
RST>? -1 

PERCE 
C H R N G 
LEN6T 

N D PRE 
.02125 
. 35573 



, 25769 
. 17836 
.00983 
. 00472 



PTS 
, 6 

NT 

E IN 
H 

V V fl L 
4 7 6 3 3 
9 1 7 1 .9 

47669 

00797 
33641 
02447 
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Converge on a Square Root 

Since most square roots are irrational, methods used to calculate them usu- 
ally involve successive approximations. Although you can simply call up the 
square root function in Basic or on many pocket calculators, it is interesting 
to explore various methods of calculating square roots without these built-in 
functions. After all, these built-in functions are nothing more than successive 
approximation routines already installed in the machine. 

Obviously, a square root is the inverse of the operation of squaring a num- 
ber. All of the methods of calculating square roots use this fact, but the way 
in which it is employed is quite different in various calculators and 
computers. 

The program here calculates an upper and lower limit for the square root 
of a number and, by successive approximations, pinches the root to within a 
smaller and smaller interval until it reaches the desired level of accuracy. 

The starting value for the lower limit is and for the upper limit the num- 
ber, Z, whose square root is sought. The program then divides this interval 
into ten steps by simply dividing the difference between the numbers by 10. 
The variable I is increased from the lower limit to the upper one by the value 
of the step, S. At any point, if I squared becomes greater than Z, a new upper 
limit is set to I and a new lower limit is set to I - S. 

This method converges very quickly and adds approximately one decimal 
place of accuracy with each pass beyond the third. What happens when you 
enter into the program a number that has an exact square root such as 25 or 
49? Why? 

Another approach to calculating square roots by successive approxima- 
tions is to start with a trial root, X. If X * X is less than the original number 
N, then increase the trial value by a 0. 1 . If X * X is greater than N, return to 
the previous value. This is the first digit of the root. Now, start advancing by 
0.01. Continuing in this way, one digit is developed at a time until the desired 
precision is reached. 

This method is quite suitable and fast for square roots of numbers less than 
1 . A good first trial root value is 0. 1 . But is it suitable for larger numbers? 
How should a starting trial value be determined? In this method, especially 
for numbers greater than 10, the initial trial value for the root matters a great 
deal in determining the length of time it will take for the calculation to con- 
verge. Write a program using this method and compare the speed and ac- 
curacy with the program in the book. 
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NT "CflLCULftTE SQUARE R 

(I IJ 

NT "SET fl C C U R n C Y lU LI 
NUMBER = " : Z 

LIMIT", " 



U T " Y n U R 

1 

NT " ":■' LOWER 

LIMIT" 



B- 

NT 
R 
R 

XT 

E* 

TO 

I 

I - 

TO 

I N 

HE 

IN 



R > .' 1 

R •> E 
B S < fl * B - Z >< E THEN 
I=R TD B STEP S 
<:i*I THEN 160 

I 
liD 

SO 



190 



•8 



" " I "RCCURflv:;Y L IMI 
"RVERRGE VALUE ="; 



U L R T E SQUARE ROOTS 

fi C C U R fl C Y IN LINE 40 
NUMBER =3 



LOWER LIMIT 


1 . 5 
1.71 
1.731 
1.7319 



UPPER LIMIT 
3 

1 . 8 
I . 74 
1.734 
1.7322 



RC 
RV 



CU 
ER 



R R C Y LIMIT PERCHED 
RGE VRLUE = 1.73205 



** DONE ** 



>RUM 
CRLCULRTE 



SQUARE ROOTS 



SET ACCURACY IN LINE 4 
YOUR NUMBER =54 



L U E R 




LIMIT 



5 



7 

7 
7 



4 

2 

344 

344 

3 4 3 :]; 2 

343428 

3 4 8 4 6 5 8 

3 4846904 



UPPER 
54 

10.8 
7.56 



LIMIT 



398 
3 4 94 
3 4 8 8 6 
3 4 S 4 S! 2 
34S4712 
3484 6 953 



RCCURRCY LIMIT REACHED 
RVERRGE VALUE = 7.3 48469 3 1 



DONE 
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Compounding 



As with successive approximations and recursion, compounding requires 
many repetitive calculations. Some compound interest and growth situations 
can be represented by a formula, but for many problems, solving by repetitive 
calculations is an excellent approach. 

There is nothing magical about compounding; you generally start with an 
initial amount of money, number of animals, etc. This quantity then grows or 
diminishes by a certain percentage at set intervals. The new amount is the old 
amount increased or decreased by this percentage. Repeat this calculation 
over and over again, and you have solved the problem. 

Five different types of problems which involve compounding of some sort 
are described in this chapter. Try the extra problems that follow some of the 
programs; you may be surprised at the results. 
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Indians and Interest 

Here is a simple compound interest problem that produces an astonishing 
answer. The problem has to do with the sale of Manhattan Island to the 
Dutch for approximately $24 worth of trinkets and beads. 

If the $24 that the Indians received in 1626 had been deposited in a bank 
paying 5-%% interest compounded annually, how much would it amount to 
in 1983? 

The program here solves this little problem by making use of the formula 
to calculate compound interest. In particular, if P dollars are invested at an 
interest rate of R (expressed as a decimal) and compounded N times, then the 
total amount A is given by the formula: 

A = P(l + R)N 

How much was gained in 1983 alone? How much was gained in the decade 
from 1974 to 1983? You can change the value of N in Statement 30 to get the 
answers to these questions. But a better way might be to enter the ending year 
with an input statement. 

Can you read a number expressed in the E (exponential format)? In a more 
conventional format, the number is $11,176,500,000 or $11.2 billion. 

There are many other ways you can improve the program and make it 
more suitable for general purpose compound interest calculations. Modify it 
to accept any starting and ending year, any rate of interest, and any starting 
principal amount. 

The problem as stated is somewhat unrealistic since banks were not paying 
interest rates of 5-%% from 1626 to 1983. Change the program to calculate 
the total amount based on the following interest rates: 

Years Interest Rate 

1626-1830 1.5% 

1831-1870 2.0% 

1871-1910 3.0% 

1911-1921 3.5% 

1922-1929 6.5% 

1930-1940 2.3% 

1941-1945 3.5% 

1946-1960 5.3% 

1961-1980 6.5% 

1981-1983 9.5% 
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P = 2 4 

R = . 5 7 
N = 1 9 8 3 
fl = P * < < 
PRINT " 
I N D I fi H S 
FDR MflfJ- 



RD BEEN 
60 PRINT 

BY 1 9 !5: 3 , 
n U N T E D T D 



-162i6 
1 + R ) '•• M > 

"IF THE S24 WHICH T 

~ ! " R E C E I V E D IN 16 2 

:"HflTTRN ISLAND H 

N V E S T E D R T 5.75-/., 
"IT b,l D U L D HAVE R M 
" S " 1 R 



■ R U N 
IF THE $24 WHICH THE INDIRMS 
R E C " ' 
HRT' 
I tJ V I 
IT ■ 
? 1 



:EIVED IN 1 
TfiN ISLAND 

'ESTED AT 5 
l^m U L D HAVE 
1 1 76 5E+10 



z t. FDR M fl N - 
HAD BEEN 
7 5 ;'. r BY 19 8: 
R MOUNTED TD 



* * Ii □ N E * * 
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Systematic Savings 

In the previous program, compound interest was calculated by means of a 
formula well known to bankers, money lenders and real estate agents. How- 
ever, if you did not know the formula, how would you calculate interest on 
money in a savings account by hand or with a calculator? 

You would probably begin by multiplying the principal amount P by the 
interest rate R and adding that amount to the original principal at the start of 
the second year. Doing that for all the years the money is in the bank will 
yield a final amount. Why not write a program to perform the calculations in 
this manner rather than use a formula? Here is such a program. 

It is set up for an initial principal amount of 100 (Statement 50), an interest 
rate of 10% (R = 0.1 in Statement 60) and ten years (N = 10 in Statement 
70). Naturally these values could be read in using input statements. The 
clever calculation in Statement 100 rounds off the amount to two decimal 
places (dollars and cents). 

In contrast with the program in the previous section, this one does not use 
a compound interest formula, but simply adds the interest each year to the 
growing principal amount in Statement 90. 



10 PRINT "CRLCULFITES INTERES 
T D N S 1 C " : " I N V E S T E r R T 1 J: 
FDR 1 Q Y R S " : " " 

2 PRINT " fi T E N D " , '" C U R R E N T " 

3 PRINT "DF YERR" . " BflLflNCE" 
40 PRINT 

5 P = 1 CI 

60 R= . 1 

7 N = 1 

80 FDR 1=1 TD M 

90 P=P+P*R 

100 PRItJT I,INT<;P*100+.5>.'liD 



110 NEXT I 



■ R U N 

CRLCULRTES INTEREST DH SlOO 

INVESTED R T 1 :■: FOR 1 Ci Y R S 

RT END CURRENT 

DF 'T'EflR BRLRNCE 

1 110 

2 121 

3 133.1 

4 146.41 

5 161.05 

6 177.16 

7 194.87 

8 214.36 

9 235. 79 

10 2 5 9.37 
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Now, how could this program be modified to allow for a plan of systematic 
saving? In other words, instead of letting the $100 lie around all lonely while 
it is compounding, each year you add another $100 to it. With this program, 
making the modification for systematic savings is easy: Statement 105 is 
added to add the new deposit each year to the ever growing principal. 

10 PRINT " C R L C U L a T E S I N T E R E S 
T D N S 1 " : " I H V' E S T E n AT 1 'v: 
E V E R Y V E R R " : " F" D R 1 V ERRS': " 

n 

2 P R I N T " R T E N D C 
U R R E N T " 

3 PRINT "DF VERR INVESTED E 
R L R N C E " 

4 PRINT 

5 P = 1 
5 5 C = 1 Ci 
60 R=. 1 

7 N = 1 

80 FDR 1=1 TD M 
90 F' = P+P*R 

100 PRUJT I;TRB(;il>;I*C;TRB( 
2 0); U4 T t P * 1 Q + . 5 .> --■ 1 
105 P=P+C 
.110 NEXT I 

> RUN 
CRLCULRTES INTEREST DN SIGO 
INVESTED AT 1 J; EVERY YERR 
FDR 10 YERRS 



RT END 




C U R R E t^ T 


D F Y ERR 


INVESTED 


B ft L R N C E 


1 


100 


1 10 


2 


2 00 


231 


3 


3 00 


364. 1 


i 


4 00 


5 10.51 


5 


5 00 


6 7 1.56 


g 


600 


843. 72 


7 


7 00 


4 3.59 


9 


8 00 


: 257. 95 


9 


9 00 


: 493 . 74 


1 


10 


17 5 3. 12 



** DOME ** 
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In an effort to attract depositors, many banks over the last 20 years have 
started offering savings and investment accounts in which the deposits are 
compounded more often than annually. In the early 50's, many banks went to 
quarterly compounding; in the late 50's, to daily; and in the early 60's, some 
"competitive" S&L's went to continuous compounding. 

How much does more frequent compounding really mean? Modify either 
of the two programs to compute the interest for more frequent compounding. 
What is the difference in interest for one year for annual, quarterly, monthly, 
daily, and continuous (every second) compounding on a principal amount of 
$1000 invested at 8%? How about for ten years? 

Incidentally, if you want to use the formula method in the previous section 
for this calculation, the formula for P principal invested at R rate com- 
pounded N times per year is: 

A = P(l + R/N)N 

Try this problem which uses the same principles of compounding. Con- 
sumer prices rose an average of 8.8% in 1980. While the government keeps 
trying to bring inflation under control, they don't seem to be meeting with 
much success. Assuming that prices continue to go up this much (8.8%) ev- 
ery year, how much will a $6000 economy car (1980 dollars) cost in the year 
2000? How much will it cost when you are 65 years old? 
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Systematic Savings Revisited 

Using the formula for compound interest and systematic savings, this pro- 
gram will calculate the amount accumulated after a given period of time. 

When you run the program, it will ask how much you wish to save each 
month, the number of compounding periods in a year, the interest rate, and 
the length of time you wish to continue your systematic savings program. The 
program will then calculate the total amount at the end of that period. 

From the preceding programs, you should be able to see that a systematic 
savings program is a very effective way to accumulate a nest egg for the 
future. 

IQ PRINT "CnLCULRTES IN TERES 

T n N II E R L n f-J C E FDR fl SYSTEM 

R T I C " : "SAVINGS P R D S R fl M " : " " 
20 INPUT "tlDNTHLY LEPDSIT =" 
■ R 
30 INPUT "CailPDUNDINS PERIOD 

C; y* Y p = " • g 

4 INPUT "INTEREST RATE (X>!. 

X ;: > = " : R 

5 R = R .' ■ 1 

60 INPUT "NUMBER QF YERRS =" 
: N 

7 PRINT 

8 P R I N T " R T E N D C 
U R R E N T " 

SO PRINT "DP YERR INVESTED B 

R L R N C E " 

10 P R I r^j T 

1 10 C=12*R 

12 P = 1 2 * R 

130 FDR 1=1 TD N 

14 P = P * < 1 + R / E ;"^ B 

150 PRINT I;TRBc: 11)51* C;TRB( 

2 ;■ ; I N T C P * 1 + . 5 > .' 1 

16 F = F' + C 

17 UE>\J I 



C R L il; U L R T E :=: INTEREST R N D 
ERLRNCE FDR R SYSTEMATIC 
SRVINC5S PRDGRRM 

MONTHLY DEPOSIT =10 
C D M P n U N D I ^J 6 P E R I D D S /- Y R =12 
INTEREST RATE <■'. K )<.>'■ >: ) =15 
NUMBER DP YERRS =2 

RTEND CURRENT 

D F YEAR INVESTED B R L R N C E 

1 120 139.29 

2 24 300.9 7 

** DDME ** 
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CFILCULflTES INTEREST AND 
BRLRHCE FDR ft SYSTEMATIC 
S R V I N G S P R D 6 R fl M 

MDNTHLV DEPaSilT =10 
CDriPDUNDING PER I DBS'- YR =1 
n«l T E R E S: T ROTE C X K . X J: > =3.5 
NUMBER OF YERRS =8 

fl T E N D CURRENT 

DF YERR INVESTED BflLRNCE 

1 1 20 130.2 

2 240 271.47 

3 3 60 4 2 4.74 

4 4 90 591.0 4 

5 6 00 771,48 

6 720 9 6. 7. 26 

7 840 1179.6 3 

8 9 60 1410.15 

* * n D N E * * 

CRLCLILRTES INTEREST AND 
ERLRNCE FDR R SYSTEMATIC 
SAVINGS PRDGRl=lM 

MONTHLY DEPOSIT =10 
C D M P G U H D I N G P E R I D D S ^ Y R =365 
INTEREST R R T E < X X . X K> = S . 5 
NUMBER DF YEARS =8 

ATEND CURRENT 

DF YEAR INVESTED BALANCE 

1 12 13 0.64 

2 240 272. IB 8 

3 3 60 4 2 7.73 

4 4 S; 5 9 6.32 

5 6 00 779.86 

6 7 20 9 7 9.68 

7 840 1197.23 
e 9 6.0 1434.08 

** DDNE ** 

Try to combine the things that you have learned from the programs in 
these sections to write a completely generalized program for systematic sav- 
ings. It should accept the following information as input: 

• Initial deposit 

• Frequency of periodic deposits 

• Amount of each deposit 

• Interest rate 

• Interest compounding frequency 

• Length of time of savings program 

Your program should produce output in tabular form showing years (or 
other time periods), amount invested without interest, total amount with in- 
terest, and the interest alone. 
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Loan Payments 

Although saving money in a systematic way is a noble goal, many people 
frequently find themselves talking to a different bank officer, namely the one 
in charge of loans. 

This program will calculate the payments for a loan for a period of one 
year or longer. The program asks you to enter the key facets of the loan in 
question: amount borrowed, annual rate of interest, interval between pay- 
ments and term of the loan in years. 

The sample run shows a relatively short-term loan (2 years) of $3000 at a 
bargain basement interest rate of 8%. The monthly payment of the loan is 
computed to be $135.68 and the total interest $256.34. Why is the total in- 
terest not $240 (8% of $3000)? Instead it is 8.54% Is this a mistake in the 
program? 

Try the program for some longer term loans at real world interest rates. 
For example, run it for an automobile loan of $8000 at 12% over a 5-year 
period. Perhaps you can see from this that systematically saving your money 
and paying cash for an item makes more sense than time payments. 

Run the program to calculate the mortgage payments on a $ 1 50,000 house 
with a $40,000 down payment. Try it with a 16% interest rate stretched over 
a 30-year period. Ouch! Look at all that interest! 

At the heart of this program are Statements 150-180. Can you see what is 
happening in these calculations? 



■ R U U 

PRDGRflM CRLCfJLRTES PRVtlENTS 
□H R LONG-TERM LORN 

fl M Q U N T E D R P. a bJ E D =12 5 
INTEREST RATE C X X . X J; ) =17 
PRVriEHT INTERVRL <riDS.> =1 
TERM DF LDRN <^'ERR'5) =25 
OUTPUT RS TREiLE <EHTER 1> 
DR TQTRLS DNLV ENTER 2>? 2 

TDTRL iriTERE:=;T = 414153.^. 1 
P R I N CI P R L R E P R I D = 12 5 
MDMTHLY PRYMErJT =S 1797.18 
RND TDTRLS S 53915.3.61 

* * Ii D N E * * 
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■ R U H 

PRDGRfltl CflLCULflTES PflYMErrrS 
DN fl LDNG-TERH LDRN 

RflDUNT EDPRQWED =3000 
INTEREST R R T E <.'K)i.. X ■■•; > = 8 
F' fi V M E N T I N T E R V R L Ol D S . > =1 
TERN OF LDRN ■: Y E fl R S ) =2 
OUTPUT R S T R E L E 'CENTER 1 ) 
DP TDTRLS DMLV EUTER 2:>? 2 

TDTRL INTEREREST = 25e. .34 
PRINCIPRL REPRID = 3000 
MONTHLV PRYMENT =5 135.68 
RHD TDTRLS S 3256.3 4 



♦* DOME ** 

PRDGRRfl CRLCULRTES PR' 
DN R LONG-TERM LORN 



' n E N T S 



R ri Q LI H T B D R R a W E D =8000 
INTEREST R R T E t! X X . X' •'. > =13 
PRYMENT INTER VRL <MDS.) =3 
TERM D F LDRN C N' E R R S ) =5 
OUTPUT RS TRELE <ENTER 1> 
OR TOTRLS ONLY ENTER 2>? 1 



T R B L E IN 3 M D N T H I N T E R V R L S 



PER- 
IOD 



P R N C P L 
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INT 
END 



DUE 
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P R N C: F' L 
REPRID 
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1 
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. 1 3 

. 91 



3 4.64 

8 13 

MONTHLY PRYMENT =S 550.23 
RND TDTRLS I 11004.64 

** BDNE ** 
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Interest on Credit Purchases 

AH too frequently, the rate of interest to be charged on a loan or credit pur- 
chase is hidden in the small type. After all, the bank or car dealer or finance 
company wants to convince you that you can afford that new car or home 
improvement or whatever. 

This program calculates the interest rate on a loan given the principal value 
of the loan, the number of payments, and the amount per payment. To make 
things even easier, the program will accept the cash purchase price of the arti- 
cle and the down payment and automatically compute the principal value of 
the loan. 

The sample run shows a rather unrealistic loan for an item costing $88.99. 
The down payment was $10.00 and the loan is over a period of 18 months; 
each monthly payment is $4.85. The program run shows that the actual in- 
terest rate is a modest 5.69%. 

Now try the program with a more realistic loan. A termite company in 
New Jersey advertises a total treatment for your house for only S200; just 
$29.95 down and 24 monthly payments of $1 1.95 per month. Is this the bar- 
gain that it seems? What is the annual rate of interest? 
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: " R fi T E D N CREDIT P 
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-RUN 

P R D G R R f1 
RATE ON 



CRLCULRTES 
CREDIT PURi 



I NTEREST 
H fl S E S 



PRIi:E DF ITEM =88.9 9 
DOWN P R V N E t^ T =10 
NUMBER DF PRVMENTS =1S 
PAYMENTS PER MONTH =1 
R MOUNT PER PRYMEr^T =4.85 

INTEREST RATE IS 5.69 'i 

*# DONE ** 
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Population Growth 

So far, all the compounding problems in this chapter have involved 
money — interest, savings, and loans. But many fascinating compounding 
problems do not involve money. Consider the following problem. 

In 1960, the population figures for the United States and Mexico were 180 
million and 85 million respectively. The annual population growth rate for 
the United States was 1.23% and for Mexico 2.23%. If these growth rates 
remain steady, in some distant year the population of Mexico will exceed that 
of the United States. In what year will this occur? 

The program uses the method of accumulating the principal amount rather 
than a formula, and applies it to both populations. The sample run reveals 
that the population of Mexico will overtake that of the U.S. in the not too 
distant future. 

10 U = 1 S 

2 M = 3 5 13 O 

3 R 1 = 1 . I 2 3 
40 R2=l . 0223 
50 Y=1960 

6 U = U * R 1 

70 M=M*R2 

S Y = Y + 1 

90 IF M<U THEN 6 0, „ ^^^ 

10 PRINT " Y E R R = " ; Y ! " POP 

ULflT r[]N" : " " , ,, 

110 PRINT " U . S . Fl . " » U 

12 PRINT " M E X I C D " , t1 

'yERR = 2037 PDPULRTIQN 

U . S . fi . 46 1 406536 . 3 

MEXICO 4644646 9 3.9 



But what if the Mexican people were diligently trying to bring their 
increasing population under control and their annual growth rate was 
increasing by only 0.001% per year, compared to the increase of the U.S. 
growth rate of 0.01% per year. If that were the case, would the population of 
Mexico ever exceed that of the U.S.? 

Insert Statements 60, 70, and 85 into the first program, run it and find out 
the answer to the question. 



>60 IJ = U*<:i.01*Rl) 

>70 M = M*';i.001*R2:» ,„ 

>S5 IF Y>9999 THEN 100 
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Going back to the first program, run it with the population data from the 
U.S. (180 million) and California (15.7 million) from 1960. At that time, the 
annual growth rates were 1.23% and 3.7% respectively. Running the pro- 
gram will indicate the year that the population of California will exceed that 
of the United States. This, of course, is nonsense. Where is the discrepancy? 

This program might make more sense if it were restated to ask in what year 
the population of California would exceed that of the remainder of the U.S., if 
ever? 

Compounding can also be used to solve other kinds of population growth 
problems. Consider the bristleworm. The bristleworm can reproduce by split- 
ting itself into 24 segments, each of which grows a new head and tail. What is 
the maximum number of bristleworms that could be obtained in this fashion, 
starting with only one worm, after ten splittings? Assuming a splitting occurs 
every 22 days, how many offspring will one worm produce in a year? When 
will the earth be overrun with bristleworms? 

Here is another problem for which the principle of compounding is useful. 
It takes nature about 500 years to produce one inch of topsoil. Many years 
ago, the United States had an average depth of almost nine inches of this 
good dirt, but as of 1975, the country was down to about six inches. This type 
of soil is necessary, of course, for growing food. 

Careless management of our soil causes about 1 % per year to erode away; 
it is then lost forever. Once soil depth reaches three inches or less, it is impos- 
sible to grow crops on a large scale. Write a program to calculate the year in 
which the U.S. will have less the 3" of topsoil, assuming that it continues to 
erode away at 1 % per year. 

Will you be alive then? Will your children be alive? Will anyone be alive? 
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Probability 



Statistics and probability are subjects that bring up all kinds of images. 
Some people who have not had a pleasant time in a required college statistics 
course keep as far away from the subject as possible. To other people, statis- 
tics is something with which devious manufacturers can mislead you about 
their products. 

Not long ago, one foreign automobile manufacturer boasted that 90% of 
their cars sold in the United States in the last seven years were still on the 
road. This sounds like excellent reliability. But consider the fact that this 
manufacturer was rapidly expanding in the U.S. market and 65% of the cars 
they had sold in the U.S. had been sold in just the previous two years. You 
would expect that all of these would still be running. 

If the manufacturer had sold 10% of their 7-year volume in Years 1 to 3, 
and most of these cars were not running, then the advertising claim loses 
much of its meaning. 

When approached in a logical, step-by-step manner, statistics and probabil- 
ity are not difficult. In fact, they can be a great deal of fun. 

Some of the programs use a formula while others simulate the event in 
which we are interested. The results are the same, but you may find that the 
simulations help you to understand exactly what is happening. 
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Pascal's Triangle — Calculated 

Pascal's Triangle is quite fascinating. As you can see from the portion of 
one reproduced below, each row is symmetrical. Each row also contains the 
coefficients for a binomial expansion. The sums across the ascending diag- 
onals form the Fibonacci sequence. The sums across the rows are all powers 
of two. Each row corresponds to the digits of a power of 1 1 . Every element is 
the sum of the two above it. And, in case you care, all the elements in it are 
identities in combinatorial theory. 

The ways this marvelous triangle can be generated are as varied and in- 
teresting as its properties, though perhaps more difficult to figure out. Here is 
one way to use a computer to generate the triangle. 

Any element can be found be adding together the two elements immedi- 
ately above it. The program uses this principle to produce a triangle eight lev- 
els deep. Lines 10-30 set the rightmost diagonal to 1. Each element is stored 
in the two-dimensional variable P(R,C) with R denoting the row and C 
denoting the "column." The variable T (Line 50) simply leaves some blank 
spaces so the output resembles a triangle. The crux of the calculation is in 
Line 70 in which the value of each new element is calculated. 

There is another interesting way to calculate Pascal's Triangle in even 
fewer statements than the program here. It generates the triangle one element 
at a time and does not use any arrays or two-dimensional variables. Can you 
figure out how to write such a program? 
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>RUN 



FDR C=l TO 9 

P < C 7 C> = 1 

NEXT C 

FDR R=l TD 9 

T=l?-R*2 

FDR C=2 TD 

P C R + 1 ■ ■ - 



\ — i. I — r-.-' ^ 
FDR C=2 TD R+1 
P C R + 1 , C .^ = P < R , C 3 + P C R , C - 1 > 
IF P ■; R , CO = THEN 110 
PRINT TfiB(:T> ; P<R, C) ; 
T = T + 4 

hi C- tJ T I-- 



NEXT C 

p R I r^ T 

P R I t>l T 
NEXT R 



1 1 

1 2 1 

13 3 1 

14 6 4 1 

1 5 10 10 5 1 

1 6 15 2 15 6 1 

^7 21 35 3 5 21 7 1 

* * B H E * * 
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Pascal's Triangle — by Probability 

This program simulates the dropping of balls through a triangular array 
shown below. At each level, a ball is equally likely to fall either to the left or 
right. At the bottom of the array a cup is placed at each end point; these cups 
collect the balls. After each group of balls has been dropped, the number of 
balls in each cup is tallied and displayed. 
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> R U N 
PRSCflL'S TRIANGLE BY 
PRDBflBIL ITY 

BALLS TD DROP =1000 
NUMBER DF LEVELS =4 

PDSITIQN BALLS 

1 51 

2 229 

3 40 7 



i 



240 
7 3 



** DONE ** 



The number of balls landing in the various cups at each level, when re- 
duced to the lowest common denominator, should approximate the numbers 
in Pascal's Triangle. Do you know why they should? 

Try a few runs of the program with a different number of balls. Do the 
numbers obtained really approximate those in Pascal's Triangle? If, at each 
dividing point in the array, every other ball went in the opposite direction, 
then Pascal's Triangle would be produced. Try running the program with 4 
balls at Level 2; with 8 balls at Level 3; and with 16 balls at Level 4. Do your 
results look like those of the previous program which calculated the triangle? 
Does the approximation come closer as the number of balls is increased? 

How can you determine how close your results from this program are to 
the exact value of Pascal's Triangle? One way is to take the number of balls 
that dropped into each cup on a given level and divide that by the total num- 
ber of balls divided by the theoretical sum of the row. So, in the sample run at 
Level 4, you would divide the balls in each cup by 62.5 (1000/16). You then 
compare that to the "correct" number and compute the percent difference. 

Here is the result of this procedure for Level 4 (actually the fifth row) in 
the sample run: 



Cup 


Balls 


^62.5 


Correct 


Deviation 


1 


51 


0.82 


1 


18.00% 


2 


229 


3.66 


4 


8.50% 


3 


407 


6.51 


6 


8.50% 


4 


240 


3.84 


4 


4.00% 


5 


73 


1.17 


1 


16.80% 



Why are the outside cups "off' by a greater percentage than those closer to 
the center? Does this also happen on other levels with different numbers of 
balls? 
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Common Birthdays 

Here is an interesting little problem in probability. In a group of ten people 
selected at random, what is the probability that any of them will share the 
same birthday? How about a group of 20 people? Of 50 people? 

Conversely, how many people would you need in a group such that there is 
a 50% probability that at least two of them have the same birthday? How 
many people would be needed for a 90% probability of an overlap? 

Try to answer these questions, either by guess or by calculation, before you 
look at the output from the program. 

This program provides a painless introduction to the world of statistics. 
The calculation is actually quite trivial. The probability that any person in a 
group has a birthday on January 1 is 1/365. If our group has only two people 
in it, the probability that both of them have a birthday on January 1 is 1/365 
times 1/365, or a very small number indeed. The probability that they have a 
common birthday is 365 times the very small number just obtained, or about 
0.27%. 

However, if there are three people in the group, the probability goes up 
slightly. Call the three people Betsy, Ken, and Larry. From the reasoning 
above, we know that there is a 0.27% probability that Betsy and Ken have 
the same birthday; also a 0.27 % probabihty that Betsy and Larry share a 
birthday; and finally, a 0.27% chance that Larry and Ken have a common 
birthday. Hence, the total probability for a group of three people is three 
times the probability for just two people. 

A group of four increases the probability over that of two people by a fac- 
tor of six, five people by a factor of 10, six people by 15, seven people by 21, 
and so on. There is a progression here, but the probability can also be cal- 
culated by the formula: 

P = 1 365-N 



365 



Can you see why this formula produces the same result as the description 
in words and associated progression above? 

Consider all the presidents of the United States. (How many have there 
been to date?) Two of them, James Polk and Warren Harding, were born on 
November 2. Is this to be expected given the size of the group? 

Since birthdays can be predicted, at least statistically, it ought to be pos- 
sible to predict deaths as well. It is interesting too that John Adams, James 
Monroe, and Thomas Jefferson all died on July 4th. Millard Fillmore and 
William Taft both died on March 8. What is the probability of that set of 
events? 
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Now it is your turn to write a program. Here is a game to be played among 
your friends. You make a bet that they have similar preferences in colors. 
Each person in the group puts up one penny and you, because you are so sure 
of yourself (and probability), put up an amount in cents equal to the number 
of people in the group. If you win, you get to keep all of the money; if your 
friends win, each of them gets double his original bet, or two cents. 

Then each person selects a color (from a larger list than there are people, of 
course). If any two have picked the same color, you win; if all have picked 
different colors, they win. 

Since you would like to win, you need to know how many colors should be 
on the Hst for different size groups to give you a better than 50-50 chance of 
winning. You might like the probability to be up around 70% or so. You can 
produce the necessary table with a five-line program. Go to it! 
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Coins in a Pocket 

The next two programs were originally written by Glenda Lappan and 
M.J. Winter and appeared in Creative Computing magazine. They are marvel- 
ous simulations for illustrating various aspects of probability. 

Coins in a Pocket is a simulation of the following situation. A newspaper 
costs 5 cents. A customer has 5 pennies and a dime in his pocket and offers to 
pay for his paper by letting you, the vendor, select at random two of the six 
coins. If you and this customer repeat this procedure for the 20 working days 
of a month, how much more or less than $ 1 .00 (20 days x 5 cents) are you 
likely to have collected? 

The program below solves this little problem. The first random coin is se- 
lected from a group of six (Line 70 in which X = INT(RND*6)). The value 
of X can be 0, 1, 2, 3, 4, or 5. If it is 0, we assume the dime was selected and a 
second pick is not made, because it will surely be a penny. Thus 1 1 cents is 
added to the running total in Line 130. 

If the first coin is a penny (X = 1, 2, 3, 4, or 5) then we make a second pick 
from the five remaining coins. Again, if the dime is chosen (Y = 0), 11 cents 
is added to the total; otherwise 2 cents is added. 

As you can see from the sample runs, after a great number of trials, the 
average amount of money collected each day seems to be close to 5 cents. If 
you would like to convince yourself that the answer is, in fact, 5 cents, con- 
sider the following. Assign a letter to each coin, A - F. Let A be the dime and 
B through F be the pennies. Since all combinations are equally probable, here 
are all the possible combinations: 

AB 

ACBC 
AD BD CD 
AE BE CE DE 
AF BF CF DF EF 

There are five combinations with a dime (5x11 cents) and ten combina- 
tions of only pennies (10 x 2 cents). Add it up and you have 55 cents plus 20 
cents divided by a total of 1 5 combinations which equals an average value of 
75/15 or 5 cents. 
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1 ID P. fl N D D f1 I 2 E 

2 ID PRINT " P R D G R fl M :=. I M U L fl T E S 
TfiKING Tl,)0":"CaiNS AT RANDOM 

F R n M fl " : " D I M E AND 5 P E f^ N I E S 



PRINT 
INPUT 



"NUMBER DP TRIALS = 



U = 

FDR K = 1 

;< = I N T < R 

IF X = 



TO N 
ND*6 
THEN 



Y=INT<RND*6> 



130 



IF Y=G THEN 130 

U = U + 2 

G a T D 14 

U = U+ 1 1 

NEXT K 

PRINT 

PRir-JT "FlVERfl6E VALUE DF 



INS IS 



U.'-N 



; E N T S " 



PRDGRflM SiriULflTES TAKING Tl?D 
C: D I N S AT RANDOM FROM fl 
DIME AND 5 PENNIES 



NUMBER D F 



TRIALS =100 
D F COINS 



AVERAGE VALUE 
4.93 4 CENTS 

* * DONE * * 

> R U M 
PROGRAM SIMULATES TAKING TWO 
COINS AT RANDOM FROM A 
DIME AND 5 PENNIES 



NUMBER DF TRIALS =100 

AVERAGE VALUE OF COINS 
4 . i3 8 CENTS 

** DONE ** 



I ^; 
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Baseball Cards 

In statistics and probability, it is frequently necessary to predict the num- 
ber of trials on average until one is successful. For example, if you are trying 
to roll a four with one die, on average how many tries will it take until you do 
SO? 

Since a die has six sides, the probability of rolling any number between 1 
and 6 is p = Ve. Thus, on any given roll, you have a Ve chance of rolling a 
four and Ve chance of rolling something else. This failure is designated q. So 
to be successful in rolling a four, we have a '/e chance on the first roll. On two 
rolls, the probability is two trials x the probability of one failure, then success 
(2 X Vb xVe). On the third roll, the probability of success is three trials x the 
probability of two failures followed by one success (3 x Ve x Ve x Ye). 

Continuing this reasoning leads to the formula for the expected number of 
trials to success: 

E=lp + 2qxp + 3qxp + 4qxp+ ... 

This series can be reduced and solved for E which leads to: 

E= _L = -L 

l-q P 

So now the answer to the original problem can be calculated. The expected 
number of trials to roll a four is 1/p = 6 tries until success. 

But there is another way to approach this type of problem with the assis- 
tance of the computer. Consider the problem. Assume there are ten different 
prizes in Crunchies cereal boxes. How many boxes of cereal would you have 
to buy to obtain the complete set? The formula above can be expanded to 
solve this problem: 

N, N ,N, ,N_^.,,1, 1,...,1 

For a value of ten, you can solve this formula by hand. However, let's say 
you would like to solve this problem for baseball cards as found in packs of 
bubble gum. If there are 50 cards in a complete set, how many packs of gum 
must you buy, on average, to get a complete set. Write a computer program 
using the general formula above to solve for a set of any size. The answer for 
50 cards is 224.96 packs of gum; how many would you have to buy for a set of 
100 cards? 

There is another way of approaching this problem. This method is similar 
to that used to randomly select coins out of the pocket. In this case, the ran- 
dom group is set equal to the total number of cards; this value is accepted as 
input in Statement 20. Each purchase is set equal to a random number be- 
tween 1 and N (Statement 150). This card is then put into its proper place in 
the collection (Statement 160). However, if there is already a card there from 
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a previous purchase, we simply increment the purchase counter (Statement 
190) but do not get any closer to obtaining a complete set. After each pur- 
chase, we test to see if the set is complete (Statement 200), otherwise we go on 
buying more packs of bubble gum. 

When the set is complete, the number of packs of gum are tallied up and 
printed out. After a set of trials, the average is computed. This averaging 
value should be reasonably close to the value obtained by the formula al- 
though it will take a great number of trials before the two numbers are within 
1 % of each other. 

I ID PRINT "PRDGRflll SIMULATES 

E IJ V UJ G '" : " B U E E L E GUM WITH B fl S 

EEflLL" 

2G INPUT "CARDS. CARDS IN SE 

R I E S = " : N 

3 ID IF |s|< 1 1 T H E tJ 6 

4 PRINT "SaRRY, 100 IS MRKI 
MUM" 

5 G a T D 2 

6 UJ P IJ T "NUMBER DP S I M U L fi T I 

U S = " : K 

7 DIM C O I' 

■3 PRINT " " : " TR I flL " , " PACKS " 

9 S; = 

i FOR I -~ 1 T D K 

110 n = 

120 FC]R J=l TD U 

13 c <: J ) = Ci 

14 HZy.l J 

15 ;< = I N T < R N D * N ) + 1 

1 fe C •:: K > = C C>? > + 1 

1 ? Q IF €■<:><> =1 THEN 190 
ISO G n T D 150 

19 D = D + 1 

20 IF D = N THEh) 220 

2 10 G a T D 15 
2 2 7=0 

2 3 FDR J = 1 T D N 

240 T=T + C<.J> 

2 5 N E y. T .J 

2 6 P R UJ T I , T 

2 7 S = S + T 

2 S N E X T I 

■ 9 PRINT " A V E R fl G E " . S .'' K 



.^ 



> R U M 
PROGRAM SIMULATES EUYIKG 

E U E E L E GUM WITH E A S E E fl L L 

C A R D S . CARDS IN SERIES =50 

NUMBER D F S I M U L A T 1 D t>l S =5 

TRIAL PACKS 

1 2S4 

'd 246 

153 
22 2 
^' 29 

fl ■•,■' E R A G E 2 3 9 



4 



* * DONE * * 
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■RUM 

PROGRRM SinULflTES BUYING 
BUBBLE GUM WITH ERSEEflLL 
CRRDS. CARDS IN SERIES =10 
NUMBER DF SIMULflTIDNS =10 



TR I RL 


P R C K S 


1 




2 3 


2 




14 


3 




28 


4 




27 


5 




1 9 


6 




28 


7 




25 


S 




2 


9 




2 6 


10 




2 2 


AVE 


RAGE 


2 3.2 



D □ N E * * 



If you run this second program for a set of 100 cards, it will sometimes run 
for a very long time before arriving at the answer. Try it with a set of 500 
cards as are used by some real bubble gum manufacturers. You can go have 
your dinner while the program computes just the first value. Be sure to 
change the dimension statement in Line 70 to C(500). 
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System Reliability 



As more and more people in the world come to depend upon mechanical 
and electronic devices in a myriad of different ways, it is important that these 
technological devices continue to function. 

Some years ago, the military came up with a measure that could be applied 
to all kinds of systems, big and small, to measure reliability. It is called 
"mean time between failures." What this means is the length of time, on av- 
erage, between breakdowns. For a tank, this may be 100 hours, while for a 
spacecraft, the MTBF must be considerably longer than the planned mission. 

As we saw in an earlier section, it is frequently desirable to break down a 
large problem into smaller subproblems. To calculate the MTBF for a space- 
craft would be quite impossible. Instead, it is necessary to start with smaller 
systems and build up to the whole. 

Consider one of the electrical subsystems of a spacecraft. It uses five 
components as shown below, two parallel systems A and B, arranged in se- 
ries. The parallel subsystems are said to be "redundant." This is one method 
of increasing reliability since the system Will continue to work if at least one 
of the components works. 



A B 

If the manufacturer of the components stated that each one has a 60% 
probability of lasting 1000 hours, what is the chance that the entire system 
will last 1000 hours? Take a guess before reading on to the solution below. Is 
your guess greater than 60% or less? Why? 

The program below is a simulation of this system. Remember in subsystem 
A, it will continue to work if any of the three components works and in B if 
either of the two components works. However, the system will fail if all three 
of the A components or both of the B components fail. 

In Statement 40, the program is set up to make 500 trials of the system. 
Statement 60 selects a random integer between 1 and 10 for each component 
for each trial. If the value of this integer is 6 or under, the component is okay 
(60% probability of working at the end of 100 hours). If it is 7, 8, 9, or 10, 
this means the component has failed, and the program prints "Bad!" 
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i 

2 
IJ 
LI 
3 
4 
5 
6 

7 

p- 

9 
1 
1 
1 
1 
i 
1 
1 
1 
1 
1 
2 
2 
22 



PP. 
E a 

S " ! 



FO 

FD 

€.<. 

IF 

C k' 

GD 

C 

N 

I 

I 

I 

P 

13 

I 

I 

P 

15 

S 



NDDMIZE 

n^l T " P R D G R fl M S I N IJ L R T E S 

F" : "CIRCUI T FDR ICOO HD 

I) II 


R 



N 

g p 



E 
R 

I N T <: 



I 

: J 
l> = 

> = 

a 

J) 

XT 

c 

c 

r: 
IN 
TD 

C 

IN 
TD 
S + 
XT 
IN 
10 



= 1 

= 1 
IN 
J> 
1 

1 1 
= 
J 
■:: 1 
<-:2 

<3 
T 

2 
<4 
<5 
T _^ 

1 
I 

T 




T D 5 

TD 5 
T k R N D * 1 > + 1 
> 6 THEN 100 



■)=1 THEN 

) =1 THEN 

> = 1 THEN 

"BAD! " ; 
2 

>= 1 THEN 

)= 1 THEN 

" BfiD ! " ; 
20 









10 
1 



" R E L I R B I L I T Y fi B D U T 
* (, S /' 5 C > + . 5 > -' 1 ; 



9 9 9 END 



P R D 6 R R M 
C I R C IJ I T 



S I n U L FI T E S USE 

FOR 1000 H □ U R 



DF 



BAD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
ERD 
BRD 
ERD 
ERD 
ERD 
REL 



ERD 
ERD 
BRD 
BRD 
BRD 
ERD 
BRD 
BRD 
ERD 
BRD 
BRD 
BRD 
ERD 
BRD 
BAD 
ERD 
IRBI 



BRD 
BRD 
BRD 
ERD 
ERD 
BRD 
BRD 
BRD 
BAD 
BRD 
BRD 
BRD 
BRD 
BRD 
BAD 
BRD 
I T Y 



* * D D r» E * * 



BAD 
ERD 
BRD 
ERD 
ERD 
ERD 
ERD 
ERD 
BRD 
BRD 
BRD 
BRD 
BRD 
ERD 
BRD 
BRD 
REQl 



1=1 D 
RD 
RD 
RD 
RD 
_ RD 
ERD 
ERD 
ERD 
ERD 
BRD 
BRD 
BRD 
EflD 
BAD 
ERD 
.IT 7 



BRD 
BRD 
EflD 
BRD 
BAD 
BRD 
BRD 
BRD 
BRD 
ERD 
ERD 
BRD 
ERD 
BRD 
BRD 
BRD 



BR 
ER 
BR 
BR 



RD 
RD 
RD 
AD 
AD 
AD 
RD 
RD 
AD 



Statements 120-140 tests whether subsystem A works and Statements 170- 
180 checks out subsystem B. If both work, a success is counted in Statement 
210. 

Now it is your turn to use what you have learned in this chapter and write 
a program to determine the mean time between failures of the system above. 
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Geometry and Calculus 



In this chapter several of the problem solving approaches from earlier 
chapters are used to solve geometric problems. You will find that many prob- 
lems can be solved in a variety of different ways — by applying a formula, by 
trial and error, by successive approximations, and, in some cases, by common 
sense. Perhaps the most important thing to learn from these problems and 
programs is how to analyze a problem to reach the solution quickly and 
painlessly. 
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Crossed and Slipping Ladders 

Here is a simple problem of a slipping ladder. A ladder 25 feet long is 
placed so its foot is 7 feet from the base of a building. The base of the ladder 
slipped on some loose gravel so that the top is 4 feet lower than where it was 
to start. How far did the foot of the ladder slip? 




(X^) 




The diagrams show the two positions of the ladder. By the Pythagorean 
theorem we know that a^ + b^ = c^, hence, the equations needed to solve the 
problem are: 

X = ^25^ - 72 



y = v/25' - (x-4)^ 
and the amount of slippage of the base is z = y - 7. These three equations are 
put into a computer program which quickly calculates an answer of 8 feet. 

While the aritmetic in the above problem is not particularly messy, it is still 
no great joy to solve by hand. However, the computer is just as happy to do 
the problem with really messy dimensions, say a ladder length of 27.83 feet 
and a distance from the wall of 7.62 feet. 



10 F' R r N T 

2 C = 2 5 

3 E = ? 

4 X = S Q R < C 

5 ■(' = S Q R < C 

6 Z = r" - E 
70 PRINT "BASE 
EET" 



"SLIPPir-IG LADDER 



2 - E '-^ 2 ) 
2 - -; X - 4 ) 



2 ) 



SLIPPED" ; z; "F 



> R U M 
SLIPPING LADDER 
ERSE SLIPPED 8 FEET 

* * HONE * * 



100 



>20 



0=27.83 
B = 7 .62 



■RUN 

SLIPPING LADDER 
EASE SLIP F' ED 3.386131551 
FEET 

** DDHE *♦ 



Let's consider an old problem found in many classic collections. Two lad- 
ders, one 25 feet long and the other 35 feet long lean against buildings on op- 
posite sides of an alley as shown below. The point at which the ladders cross 
is 12 feet above the ground. How wide is the alley? 




By using similar triangles twice we find that 

12/a + 12/b = 1 
Then, by applying the Pythagorean theorem and reducing, we obtain: 

a" - b' = 600 
Using one of the methods described in the Problem Solving chapter, you 
can solve these two simultaneous equations. Or they can be combined into 
one equation in which z, the width of the alley is 

z = V35' - a" 
By eliminating b, the following fouth degree equation is obtained 

a* - 24a3 - 600a' 4- 14400a - 86400 = 
Again, this can be solved using one of the methods in the chapter on Prob- 
lem Solving. This is the traditional approach, but there is another. The solu- 
tion is the intersection of the curves described by the two original equations. 
By using successive approximations, say in steps of 1 for a, you could solve 
for bj and b^ in the following restated original equations 
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12a 
a- 12 



b^ = ^7600 - a^ 

When b, falls below b^, reduce the step to 0. 1 to obtain a closer approxima- 
tion. By continuing this method of successive approximations, it is possible to 
obtain a very accurate solution. 

Is there another approach? Yes, there is and it also avoids the quartic equa- 
tion. It uses the original equations in a trial and error procedure as described 
in the chapter on Sets and Repetitive Trials. See if you can write a program 
using this approach. 

Here is another classic problem for you to solve in any way that you like. 
What is the longest ladder that can be carried in a horizontal position around 
the corner made where a 12-foot wide alley meets one that is 8 feet wide. The 
diagram shows the problem. 




12 
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Distance Between Coordinate Points 

This program solves for the distance between any two points in three- 
dimensional space defined by their x, y, and z coordinates. 

The formula for the distance between two points in three-dimensional 
space is: 

d = (X, - x,)=^ + (y^ - y,r + (z. - z,)^ 

It would be desirable for the program to be able to solve any problem of 
this kind. One approach would be to use input statements to accept the point 
coordinates. Another is to use a data statement to define the points. Then 
only this one statement has to be changed for a new problem or set of 
problems. 

The program is written to calculate the distance between three sets of 
points. The points used were: 

0, 0, I and 3, 4, 5 

3.5,-4.716.2 and -0.9,3.0,4.4 
67, 36, 82 and 54, 25, 90 

This calculation is used I extensively in aerospace navigation. Can you deter- 
mine the angle or "compass heading" of the resulting flight path? Better start 
with just two dimensions. 
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10 PRINT "PR 

D I STBNCE " : " 

N T H R E E - " : " D 



a G R Fl ri C fl L C U L Ft T E S 
FETWEEN PDINTS I 
in E N S I D N R L S P R C E 



2 PRINT 

3 F' R I N T 
T R N C E " : " 

4 F: E R D 

5 L = S Q R 


" C 

M 

R, E 
■:. <R 


POINT" 
QDRD INRTES 

, C , D , E . F 

- D > -^ 2 + < E - E V^ 2 


rii3 
+ <c.- 


^. PRINT 
7 PRINT 
SO PRINT 
■5 G D T D 
100 DRTR 
1 10 DflTR 
4 . 4 

120 DRTfl 
> 


4 

g , 

6 7 


B; C 
e; f; 

0,0, 
5, -4 

, 36 .■> 


IRE 

3, 4 
. 7 , 

3 2 , 


c: 1 6 ) ; L 

, 5 
6.2,-. 

54,25, 


•? , 3 , 
9 


> 

> R U N 
PRaGRfltl CRLCLILRTES 
BETWEEN POINTS IN 
DINENSIDHRL SPRCE 


D I S T R t^l C E 
THREE- 


POINT 
C n a F' D I N fl T E S 




DISTANCE 





3 4 5 






t « 


071067 


S12 


-?9^ i"- 


7 
4 . 4 


6. 2 


9. 


4 9 3 9 


36 6 


%l i% 


82 
•3 




18 


. 81488 


772 



* DRTfl ERRDR IN 40 
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Area — by Calculation 

It is a simple matter to calculate the area of regular geometric figures by 
using the usual formulae. However, the problem becomes more difficult when 
it is necessary to calculate the area of two combined regular shapes or the 
area of irregular shapes. This program shows the method of analysis for solv- 
ing a problem of the first type, while the next program demonstrates four 
methods of dealing with irregular areas. 

The problem is to solve for the shaded area of the figure below for any 
value of the radius, R. 




The first step is to recall the formulae for calculating the area of a circle 
and square: 

A (circle) = pi * R' 
A (square) = Side^ or (2 * R)' 
The difference in area between a square and a circle inscribed within its 
borders is: 

A (difference) = A (square) - A (circle) 
and the area of one corner is the difference divided by 4. The program below 
will calculate the area for any radius. 

Now it is your turn. Write a program to calculate the difference in area 
between a circle and square in which the square is inscribed within the circle. 
Now, change your program to calculate the difference in area for a triangle, a 
hexagon and an octagon inscribed within a circle. 



10 PRINT "PROeRflM 

flREfl" : "BETWEEN Fl 

C R 1 1 E D WITHIN fi 



CftLCULFlTES 
CIRCLE INS 
SQUARE" ! " 



2 
30 
40 
5 
60 
70 

3 



INPUT "R ALIUS =":R 

= 3. 14159*R'-2 

.§ = <; 2 # R > •*■ 2 

f I = S - C 

S 1 = n ./ 4 

PRINT "TRf=lPPED fiRER ="5D 

PRINT " D N E CORNER = " ? H 1 
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> R U N 

PRQ5RflM CflLCULRTES flREfi 
BETWEEN fi CIRCLE INSCRIBED 
WITHIN fl SQUARE 

RADIUS =10 

TRflPPEB RREfl = 85.841 

ONE CDRNER = 21.46025 

** DDNE ** 

> R U N 

PRDGRflM CRLCULflTES RREfl 
BETWEEN R CIRCLE INSCRIBED 
WITHIN fi SQURRE 

FiflDIUS =3.14159 

TRAPPED RRER = S. 472152802 

ONE CDRNER = 2.1180382 

** DDNE ** 



Extend your program to calculate the difference in area for any regular fig- 
ure inscribed within a circle of radius 1.0. Set up a table or results as follows: 
Number of sides Difference in Area 

3 
4 
5 
6 
7 
8 
What can you conclude from these results? Do these areas follow in some 
sort of progression? 
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Area — by Integration 



In many cases it is necessary to determine the area of an irregular figure or 
the area under a curve where an exact formula is not available. The approach 
most commonly used is to divide up the enclosed area into small regularly 
shaped pieces and sum up the areas of all of these pieces. 

The easiest shape to use in these calculations is a rectangle. A group of 
rectangles can either be inscribed within the irregular figure or curve, or 
circumscribed around it. The first two diagrams show these two methods be- 
ing used to find the area under a curve. A third method is to use trapezoids 
which, depending upon the direction of curvature, will either be inscribed or 
circumscribed automatically. 





Inscribed 
Rectangle 



Circumscribed 
Rectangle 



Trapezoid 
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A fourth method, known as Simpson's Rule, essentially fits a series of 
parabolas between the points of the curve and calculates the average area. It 
requires that the area be divided into an even number of parallel slits. Let us 
call the total number of divisions 2m which are h distance apart. The point on 
the curve where the first line intersects is yo, the second yi, and so on until 
y2m- The area is then given by the formula: 

A= VM(yo + yzJ + 4(y, + yj + ... 
+ y2m-l) + 2(y2+ y4 + • ■ ■ yarn -2)] 

The area as calculated by Simpson's Rule converges extremely quickly 
compared to the other methods. Nevertheless, as the number of intervals in- 
creases, all the methods approach the same limit. As might be expected, the 
trapezoidal approach converges more quickly than either of the approaches 
using rectangles. Compute the average of the two methods using rectangles. 
What do you get? Does this suggest another method? 

The methods used in this program involve the calculus. And you thought 
the calculus was difficult! Now you know otherwise. 



i G F' R I H T " P R D G R R 11 C: D M P U T E S R 

RER UNDER": "R CURVE BV FOUR 

METHODS F '■ ! " H U t1 E R I C R L I N T E G 
R fl T ID M " : " " 

2 PgINT " DEFINE CURVE IN LI 

30 PRINT "INPUT START RND EN 

D V R L U E S " 

40 INPUT "FDR X CSMRLLER FIR 

S T .■) = "; fl , B 

50 PRINT 

ig^PRiyj "NO. INTERVALS" :" IN 

SCR R E C T S T R fl P E Z I D 3 " : " C I 

RCUM RECTS PRRfiBDLflS" 

7 P R I N T 

SO M = - 2 

9 £; = 

1 13 D E F F N fl < X ) = J•^■^ 3 
110 fl = M + 3 

120 FDR N = f1 TO M + 2 

13 C = S 

140 Q=0 

150 P=0 

16 D = 2 ■-■ N 

17 H=(;E-fl;)/'D 

18 FDR 1=0 T D D - 1 
190 M=fl+I*H 

2 P = P + H * F r^ fl <; K > 

2 10 Q = Q + H * F N R ■; X + H > 
2 2 NEXT I 
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230 

240 

2 50 

260 

270 

2S0 

290 

300 

310 

320 

5 3 

340 

350 

360 

. OOi: 

370 




NEXT H 

IF i)< e. 4 

IF ft B S k' 
5 1 THEN 
END 



C C - S 
1 10 



>RUN 
PROGRflM CanpUTES AREA UNDER 
R CURVE BY FOUR METHODS OF 
NUMERICAL INTEGRfiTIDN 

DEFIWE CURVE IN LINE IGO 
INPUT STRRT RND END VALUES 
FDR X <;SI1RLLER F 1 R S T ) = 1 , 1 



NO. IHTERVfiLB 
INSCR RECTS 
C I R C- 1) H RECTS 



2 

75 3. 

524S 

4 

150 1 

37 48 

8 

19 6"9 

30 93 

16 

2 226 

27S8 

32 

2 36 1 

2642 

64 

2429 

2570 



1 875 
.6875 

. 17 13 7 5 
.921375 

. 13 6 7 19 
.011719 

.612305 
. 5 4 9 3 5 

. 2 2 3 3 8 9 
. 192 139 

.997253 
.481628 



TRftPEZDIDS 
PRRRBDLflS 



3 000. 9375 
2 4 9 9. 75 

2625. 046S75 
2 4 9 9.75 

2 5 3 1.074219 
2499. 75 

2 5 7.581055 
2499. 75 

2501 . 7077 6 4 
2499 . 75 

2 5 0.239441 
2 4 9 9. 75 



* * DONE * * 
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_s^ 



^a4 



"^m 




"According to the computer simulation, it should hit the earth in 0.00298 seconds. 
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Science 



Using techniques and approaches presented in the previous chapters, this 
chapter contains five programs in the area of science. One uses a formula to 
solve simple gas problems, two are drill exercises on the gas laws of Boyle and 
Charles, and the last two are simulations. You will see that the simulations 
draw upon many previous techniques such as progressions and repetitive 
calculations. 
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Gas Volume 

Here is a simple program to produce a not-so-simple table of values for gas 
volumes. 

The volume of a gas varies directly with the absolute temperature T (Kel- 
vin) and inversely with the pressure P. If a certain quantity of gas occupies 
500 cubic feet at a pressure of 53 pounds per square foot and an absolute tem- 
perature of 500 degrees, what volume will it occupy at 600 degrees absolute 
temperature and pressures from 100 to 1000 pounds per square foot in in- 
crements of 50 pounds? 

The original conditions are used to solve for the constant K in Statement 
40 (K = V*P/T). Then new volumes are computed for varying pressures 
with T equal to 600 degrees. The formula used is V = K*T/P. 

In the second part of the program, Lines 70 to 100 are replaced to produce 
a plot of the gas volume for the various pressures. 

How would you modify the program to deal with a more general case (i.e., 
other gasses and different temperatures)? Second, can you write a program 
that produces a table of values for a gas at different pressures and 
temperatures? 



10 \' = 5 

2 P = 5 3 

3 T = 5 

4 K = \' * P ..' T 

5 PRINT " P R D G R a M C D M P U T E S G 

elsuRES"^'^"?-""^ different'p^ 

6 I N P ij T " T E M P E R fl T U R E f. K > = " 
70PRINT " ": "PRESSURE", "VOL 
IP f^5f;. ?7l.§<^ TO 1000 3TEP 50 

10 PRINT P ^ V 
110 NEXT P 



PRDGRRH COMPUTES QRS VOLUMES 

AT DIFFERENT PRESSURES 

TEMPERATURE OO =400 

P R E 8. S U R E VOLUME 

10 3 18 

15 2 12 

2 15 9 

2 5 12 7.2 

3 10 6 

35 ? 0.85714286 

4 ID 7 9.5 
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4 50 

5 00 

5 50 

6 00 
650 
700 

7 50 
800 
850 
9 00 
9 5 

1 G 



70. 

63. 

57. 

5 3 

48. 

45. 

4 2. 

39. 

37. 

35. 

3 1 : 



66666.61E.7 

8 1 8 1 S 1 8 2 

9 2 3 ID 7 6 9 2 
4 2 8 5 7 14 3 
4 

75 

4 117 6 4 7 1 

3 3 3 3 3 3 ;i 3 

4 7 3 6 3 4 2 1 
8 



10 



' = 500 



*.' P = 5 3 

30 T=500 

4 a K = M * P .•- T 

5Q PRINT "F'RDGRRM CDMRUTES G 

as VOLUMES" !"RT DIFFERED PR 

E 3 S U R E S " 5 " " 

60 INPUT "TEMPERATURE <.K-J =" 

! Tl 

7 PRINT "PRES-'lOO PLDT 

OF VOLUME" 

80 FDR P=100 TD 1000 STEP 50 

9 V = K * T 1 .'■ P 

10 PRINT P /' 1 5 T ft E < 4 + I N T < V -' 
I 3 .■) ) J " * " 

110 NEXT P 



■'pRDGRflM COMPUTES SfiS VOLUMES 
AT DIFFERENT PRESSURES 



TEMPERflTLIRE 

P R E S ^ 1 
1 

1 . 5 
2 

2. 5 

3 

3.5 * 

4 * 
4.5 * 

5 * 
5.5 * 

6 * 
6.5 * 

7 * 
7.5 * 

8 * 
8.5 * 

9 * 
9.5 * 

10 ^ 



<: K :j =600 

PLDT DF 

* 
* 
* 



VOL u n E 
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Charles' Law Drill 

In the previous program, gas volumes were calculated using Charles' Law 
and Boyle's Law. Here is a drill and practice program (remember Chapter 1?) 
that produces problems relating the volume and temperature of a gas. When 
pressure is constant, the volume/temperature relationships can be stated as 
follows: 

Vo = V^ 

To T| 

The program presents four problems, one to solve for each of the four vari- 
ables in the equation above. 

How can you make the program more efficient? More interesting? 



10 RFINBDriIZE 

20 PRINT "DRILL DN CHARLES' 

L R W " : " " 

3 PRINT " V D L U M E S IN M I L L I L I 
IRES TEMPERATURES IN D E G R E 

P C t^ [| J II 41 

4 V = i N T < R: N D * 5 + 5 > * 1 

5 V 1 = I rJ T ( R N D * 5 + 5 > * 1 

6 T=INT<P:ND*125+125> 

7 T1 = INT<RND*12.5+125> 

'"' fl &" f' i 1 

1 6 D isi C: G D T D 1 > 1 4 , 1 8 , 2 2 
, 260 
10 V = 

110 PRINT "SOLVE FDR V WHEN" 
12 Q 1 = V 1 * T .'■ T i 
130 GDTD 27 
140 V1=0 

15 PRINT " W H fi T IS VI GIVEN" 
160 Ql =V*T 1 /T 
170 GDTD 2 70 
ISO T=0 

190 PRINT "Cl=lLCULflTE THE VflL 
IJ E D F T " 
200 Q1=V*T1/'V1 
210 GDTD 270 
220 T1=0 

2 30 PRINT "SOLVE FDR Tl GIVE 
N" 

240 Q1=T*V1-^V 
2 50 GDTD 270 
260 STOP 

2 7 P R I N T " V VI T 

T 1 " 

2 8 PRINT V ; T fl B C 7 ;■ ; V 1 ; T fl E C 1 4 
.^ ; T ; T fl B < 2 1 > ; T 1 

290 INPUT "YDUR PiNSUER =":Q 

300 PRINT "CORRECT VRLUE ="5 

Ql 

310 PRINT 

3 20 GDTD 40 
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■RUN 
DRILL D H CHARLES 



LAW 



VDL 
TEM 

SDL 

V 


YDU 
CDR 

I..J H n 

V 

63 
■f D U 
CDR 

CRL 

V 

6 8 
YDU 
CDR 

SDL 

V 

86 
YDU 
CDR 



UMES 
P E R fl T 

VE FD 

V 

<9 

R RNS 

R E C T 

T IS 

V 

R RNS 
R EOT 

C U L fl T 
w 

GO S 

R R rj S 

R E C: T 

'v' E F D 
V 
6 
R HNS 
RECT 



IN MIL 
U R E S I 

R V bJ H 

1 

80 

W E R =6 

V R L U E 

VI 5 I V 
1 

IjJ E R =6 

V R L U E 

E THE 

1 

900 

(jj E R =1 

V R L U E 

R Tl G 

1 

300 

U E R =2 

V R L U E 



L I L I T R E S 

N DEGREES K 



EN 

T 

169 

600 

= 6 9 

EN 

T 

152 

5 00 
= 65 2 1 

\' R L U E 

T 



6 1 
= 153 



Tl 

240 

833333 



T 1 
146 

5 7 3 9 4 7 

OF T 
T 1 

2oe 

9213483 



I V E N 

T Tl 

24 8 



= 19 6.0930 2 33 



* * DONE * * 
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Boyle's Law Drill 

Boyle's Law describes the behavior of gases under ideal conditions when 
pressure and volume are varied. When the temperature is constant, Boyle 
found that: 

Po * Vo = P, * V, 

Pressure is normally measured in centimeters of mercury while volume is 
measured in millilitres. As with the drill on Charles' Law, this program 
presents four problems, one to solve for each of the four variables in the equa- 
tion above. 

Unlike many other drill and practice programs, these two do not compare 
your answer with the correct one. Instead, they leave that up to you to do. Is 
this desirable? Why or why not? If you feel it is undesirable, change the pro- 
grams so they do compare answers and calculate a score. 



10 RRMD0MI2E 

20 PRINT "DRILL DN BDYLE'S L 

RU" ! " " 

3 PRINT "VOLUMES IH MI L L I L I 
TRES" : "PRESSURES IN CM DF flE 

p c u p y " ' " " 

4 V = I N T t; R rJ D * 5 + 5 > * 1 
50 Vl = INT<;RND*50 + 50.^* 100 
&0 P=INT<RND*150+1505 

7 P I = I N T ( R N B ^v 1 5 + 1 5 > 

SO C=C+1 

90 QN C GDTQ 100,140,130,22 

,260 

10 V = 

110 PRINT "SaLS''E FOR V WHEN" 

120 Ql=Pi*Vl./P 

130 GOTO 270 

140 V1=0 

15 PRINT " W H ft T is VI GIVEN" 

160 J31=P*V/P1 

170 GQTD 270 

180 P=0 

190 PRINT "CftLCULHTE THE VflL 

U E OF T '■ 

2 00 Q1=P1*V1/'V 

2 10 GOTO 270 

22 P1=0 

230 PRINT "SOLVE FDR 11 GIVE 

N " 

2 4 Q 1 = P * V z \' 1 

250 GOTO 270 

2 6 STOP 

2 ? P R I N T " V VI P 

PI " 
2 8 P R IfJ T V ; T fl B t 7 > ; V 1 ; T fl B '.: 1 4 

:■ ; P; tre<:21 > ; pi 

2 9 INPUT "YOUR fi N S U E R = " : Q 
300 PRINT "CORRECT VALUE ="; 
Q 1 

3 10 PRINT 

3 20 GOTO 40 
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> R U H 
DRILL 

VOLUME 
P R E S S U 

SOLVE 
V 


V OUR fl 
C R R E C 

Ul H R T 1 

V 

9400 
YOUR R 
C a R R E C 

C fl L C U L 
V 
S800 

Y D U f: fl 
C R R E C 

SOLVE 

V 

6300 
YOUR R 
C: D R R E C 



DN BDYLE 'S LRM 



a IN 
RES 

FOR 

VI 

740 
N S W E 
T VR 



MIL 
IN C 



L I L 1 T R E S 

M D F MERCURY 



V W H E U 
P 

251 
350 
= 734 1 





R =7 

LUE 



PI 

249 

3 5 8 5 7 



S VI 5 1 '.y 

\' 1 



H S U E 

T VR 



:r =1 

[LUE 



EN 

P P 1 

232 2 15 

0060 

= 10 143. 2558 1 



RTE 

V 1 

6 70 
NSUE 
T VR 

FOR 

VI 

970 
NSWE 
T VR 



THE 



R =1 

LUE 

Tl G 

O 

R = 1 
LUE 



VALUE 

P 



44 

= 146 

I V E rj 

P 

265 

70 

= 172 



OF T 
PI 
193 

9431 8 1 8 



F'l 


113402 1 



* * DONE * * 
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Photoelectric Emissions 

When light of a short wavelength falls on a metal surface, electrons are 
emitted from the metal. According to the description of this phenomenon by 
Einstein, there is a maximum wavelength for every metal above which no 
electrons are emitted. This is called the critical wavelength of the metal. 

This program simulates a laboratory experiment in which a metal is placed 
in a vacuum and bombarded with soft X-rays. The number of electrons emit- 
ted is collected and measured with a microammeter. The program simulates 
three trials at each of nine wave lenghts. 

After each set of experimental data, the program asks if you would like an- 
other run at a higher light intensity. The reason for doing this is that some- 
times at low light intensities, not enough electrons have been emitted for 
meaningful measurements. 

You can increase the precision of the experiment by increasing the number 
of wavelengths at which it is run. This value can be changed in Statement 60; 
note that the variable L is divided into 1000 in order to express the wave- 
length in Angstroms. One Angstrom (A) equals 10"* centimeters or 10' 
microns. 
Here are the coefficients for several metals: 

Silver .308 

Bismuth .338 

Cadmium .318 

Lead .340 

Platinum .385 

It is a rare physics laboratory in a high school or college today that has the 
experimental apparatus to run this experiment, yet with a small computer the 
equipment can be simulated. There are many other things with which you 
would not normally be able to experiment that can be simulated with a com- 
puter. Things such as a nuclear power plant, a malaria epidemic, an urban 
mass transit system, and a bicycle factory. 

Can you write a simulation for a real world system? You will find sections 
on simulations in the books Computers in Mathematics and Computers in Sci- 
ence and Social Studies. These, along with articles in 99'er and Creative 
Computing magazines, might be of some help in writing a simulation of your 
own. 
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10 P R I N 

M E T R L " : 
T X-RflV 
2 I N P U 
R L = " : V 
:5 K=IN 
40 PRIN 
E R E S " : " 
50 PRIN 
TRIAL 

iso for' 



70 
SO 
90 
100 

1 10 
120 
130 
) ) 
140 
150 
i 6 
170 
180 
190 
T I 
200 
210 

2 2 
TO 

2 30 
240 
250 

26 

27 



M=IN 
PRIN 
FOR 
IF 
I=S 
GQT 
I=S 

N=I 
PR! 

NEX 
PR I 
NEX 
PRI 
NT EN 
INP 
IF 
PRI 
10> 
INP 
K=K 
PRI 
SDT 
STD 



T 

"B 

S" 

T 

TC 

T 

II 

T 

II m 
m 
II 

L = 

T< 

T 

.J = 

L> 

QR 

D 

QR 

NT 

NT 

T 

NT 

T 

NT 

SI 

UT 

fiS 

NT 
i> ■ 

ut 

*F 
NT 

a 

p 



"PRO 
DtlBfl 



6RlR 
RDM 



II SItlULRTES 
ENT WITH SDF 



"COEFFICIENT DF MET 



1 + 2* 
"DUT 

"UflV 
"LEN 

U 

. 4! 



RND 
PUT 

E 
GTH 

TO 

i^L> 



) 



IN M I C R D fl M P 



TRIAL 
I 



TR I flL 
2 



25 STEP 







10 0. 

M; 

1 TD 3 

V THEN 130 

< I N T >; 2 5 * R N D > > 

140 

C K * K * 1 + I N T < 3 5 * R N D 



C 1 0* 
TAB 

J 



TV" ! 

RS 
= " N " 

"BY 



I + . 
<: J* 



5 > /■ 1 

7 :> ; N ; 



: " INCREASE LI5H 
" < V □ R N > " ; 



THEN 270 
WHAT FACTOR 



i 1 



50 
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>RUN 
PRneRRM SiriULftTES 
EDtlBRRDMENT WITH ; 



METRL 

OFT X-RRYS 



CDEFFICIEhJT OF METRL = 
OUTPUT IN MICRDRMPERES 



34 



WRVE TRIRL TRIAL TRIRL 
LENGTH 1 2 3 



2380 
2500 
26 3 1 
2777 
2941 
3125 
3333 
35 7 1 
384^, 



20 . 3 
2 0.2 
20. 5 
2 0. 8 
2. 3 
4 . 1 
3 

4 . 1 
3 



20 

2 

20 

20 

1 

1 

4 

3 

4 



4 
5 
6 
4 



7 
4 



20. 3 
2 0.2 
20 

20. 4 
4. 5 



4 
4 
1 
3 



4 
6 



I N R E R :=. E LIGHT 

< V D R N > ? Y 

BY IJHRT FRCTDR 



I N T E H S I T Y 
CI T D 1 r 



10 



IJ fl V E TRIRL TRIAL TRIRL 
LENGTH 1 2 3 



2 380 


2 00 


2 00 . 1 


20 


2 500 


200. i 


2 0. 1 


20 


2&31 


2 0. 1 


200. 1 


200 


2777 


200. 1 


2 00. 1 


20 


294 1 


4. 6 


3 . 6 


2. 2 


3125 


4.4 


3. 5 


3. 9 


3 3 33 


2. 4 


3. 3 


4. 4 


3571 


3. 7 


3. 7 


2. 4 


3846 


3 


4.5 


3. 7 


INCRER 


S E LIGHT INTENSITY 


<Y DR 


N > ■? U 







♦* DONE ** 
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Mutation of Moths 

This program is a simulation of the growth of a colony of pepper moths. 
The program allows a genetic mutation to be introduced in some year be- 
tween 1 and 30. The mutation can favor either dark or light colored moths. 

The program as it appears starts with a total colony of light moths; you 
could add an initial group of dark colored moths in Statement 35 as PI. 

The sample run shows a mutation which favors dark moths occurring in 
Year 3. The mutation affects about 2% of the light moths each year and 
causes them to become dark. The program displays the number of moths of 
each color over a 30-year period. 

Obviously, because of the long time period involved it would be difficult to 
carry out this experiment in school, so the computer again is of real benefit. 



iO PRINT "PRDGRRM SIMULFITES 
MUTflTiaN":"nF fl COLDNN' DF MD 

T LJ C; " ■ »' " 

20' INPUT "RRTE DF MUTftTIDN < 

I - 1 > = " : M 

-.3 F = i O O 

40 Z=PO 

50 PRINT "IN WHAT YEAR DDES 

i: H fl N G E " ! " I N ENVIRONMENT D C C U 

60 ' INPUT X 

70 L=] 

8 D = 2 

90 PRINT "CHANGE FflVDR LIGHT 

D R D R R K " : '■ M D T H S < L D R D ) " ; 
100 INPUT ES 

110 PRINT " ":"VEftR DARK MDT 
MS LIGHT MOTHS " 
12 FDR T=l TD 30 
130 IF T>=X THEN 160 
140 P1=0 
150 GOTO 2 20 
160 IF ES<>"Ii" THEN 140 
170 Pl = INT < P 1 + . 01 *M*PO+. 5> 
130 P0=INT<:Z-Pl + .5> 
190 IF PKZ THEN 220 
200 Pl=2 
210 P0=0 

220 PRINT T5TRB<8>;P1;TRB(19 
) ; PO 
230 NEXT T 
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>RUN 
PRDGRRM, SrriULRTES MUTflTIDN 
DF ft CDLDWY DF MDTHS 

RRTE QF MUTflTIDN <l-10-5=2 
IN WHAT YEfiR DDES CHRHGE 
IN ENVIRONMENT DCCUR-? 3 
CHANGE FflVDR LIGHT DR DARK 
MPTHS (. L QR D>? D 

YEfiR DARK MDTHS 
1 


2 00 
396 
588 
776 
960 
114 1 
13 13 
1 492 
1662 
1 3 2 9 

1 992 

2 152 

2 309 
2463 
26 1 4 
2762 
2907 
3049 

3 188 
3324 
345S 
3589 
3717 
3343 

4 087 
4205 
4321 



The so-called "killer bees" that were introduced into South America in the 
mid 70's were supposed to help honey production because they were much 
more energetic than the lazy honey bees in Brazil. The idea was that they 
would mate with the existing honey bees and produce a more productive 
strain. However, they went on a rampage killing people and terrorizing the 
country. Then they started to migrate north. U.S. agricultural officials be- 
came alarmed that they would invade this country and bring death and 
destruction. 

Over the years, the killer bees have started to mate with the more docile 
South American honey bees, but only at the rate of 3% per year. Assuming 
they pose a danger to the U.S. until their numbers are reduced to 15% of 
their original quantity, how many years will it take for this to occur? The 
present migration patterns will bring them to the U.S. by 1989; will they still 
be dangerous (assume that the year of introduction was 1 974)? 
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IGHT riDTHS 


10 


10 


9 80 


9604 


9412 


9224 


9 04 


S859 




8 5 0.1 


8 339 


8171 


8008 


7848 


7691 


7537 


7386 


7238 


7093 


6951 


6812 


66 7 6 


65 4 2 


64 1 1 


6283 


6 157 


60 34 


59 13 


5795 


5 6 79 



Projectile Motion 

The path followed by a projectile is called its trajectory. The trajectory is 
affected to a large extent by air resistance, which makes an exact analysis of 
the motion extremely complex. We shall, however, neglect the effects of air 
resistance and assume the motion takes place in empty space. 

In the general case of projectile motion, the body (bullet, rocket, mortar, 
etc.) is given an initial velocity at some angle 9 above (or below) the 
horizontal. 

If Vq represents the initial velocity (muzzle velocity), the horizontal and 
vertical components are: 



ox 



VpCos9 



Voy = Vi"^ 



Since we are neglecting air resistance, the horizontal velocity component 
remains constant throughout the motion. At any time, it is: 

Vj^= Vqj^ = Vq COS0 = constant (1) 

The vertical part of the motion is one of constant downward acceleration 
due to gravity. It is the same as for a body projected straight upward with an 
initial velocity Vq sin 9. At a time "t" after the start, the vertical velocity is: 

^y = %-gt = VoSin9-gt (2) 
where "g" is the acceleration due to gravity. 
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The horizontal distance is given by: 

^ = Vt = (V°^^)t (3) 

and the vertical distance by: 

y = V^yt-l/2gt2 

= (VpSin0)t-l/2gt^ (4) 

The time for the projectile to return to its initial elevation is found from 
Equation (4) by setting y = 0. This gives 

2 V^sin e 

The horizontal distance when the projectile returns to its initial elevation is 
called the horizontal range. "R." Introducing the time to reach the point in 
Equation (3), we find: 

2 V^^sin e cos 9 
R = (6) 

Since 2sin 6 cos 6 = sin 26, Equation 6 becomes: 

V„^ sin 26 
R = ^-^ (7, 

The horizontal range is thus proportional to the square of the initial veloc- 
ity for a given angle of elevation. Since the maximum value of sin 20 is 1, the 
maximum horizontal range, Rmax is Vo^/g. But if sin 26 = I, then 26 = 90° 
and 6 = 45°. Hence the maximum horizontal range, in the absence of air 
resistance, is attained with angle of elevation of 45°. 

From the standpoint of gunnery, what one usually wishes to know is what 
the angle of elevation should be for a given muzzle velocity vq in order to hit 
a target whose position is known. Assuming the target and gun are at the 
same elevation and the target is at a distance R, Equation (7) may be solved 
for (9. 



9 = 1/2 sin" 
= 1/2 sin" 




(8) 



Provided R is less than the maximum range, this equation has two solu- 
tions for values of 6 between 0° and 90°. Either of the angles gives the same 
range. Of course, time of flight and maximum height reached are both greater 
for the high angle trajectory. 

124 



For example, say the maximum range of our gun is 10,000 yards and the 
target is at 5,900 yards: 



1/2 sin 



_i 5,900 



10,000 
= 1/2 36° 
= 18°, or 90° -18° = 72° 

Try the computer game Gunner which appears below. Use trial and error 
to try and destroy the target (see sample run). You get five chances per target. 
How many shots did it take to destroy all five targets? Did you ever fail to de- 
stroy a target in five trials? 

From the discussion above, you should realize that 45 degrees of elevation 
provides maximum range with values over or under 45° providing less range. 

The maximum range of the gun will vary between 20,000 and 60,000 yards 
and the burst radius of a shell is 100 yards. 

You can also determine the correct firing angle from sine tables, or a slide 
rule, a scientific calculator or a computer. You should be able to destroy ev- 
ery target with just one shot. What happens when the target is very close? 
Can you always use whole angles? 

Now write a computer program to accept the maximum range of your gun 
and the range to the target and then calculate the correct firing angle. You 
will have to solve two problems to write such a program: 

1 . The Basic language does not have an ARCSIN function. However, the 
following formula may help. 



sin X = tan 



(n/M^J 



2. You must convert from radians to degrees. 
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B 
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ZE 

0000 

MAX 



R R E C D M M fl N D I M 6 
HIT WITHIN 10 
KE TRRGET":: 

*RND+20000) 
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NT<R*<.1+.8*RND5> 



D 25 
INT 
TD 2 
INT 
TD 2 
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TD 2 
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II 
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= S1 + 
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60 
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P 
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S+1 
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DM ! " ! s 
INT "ENEttY GDT YDU FIR 
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G CDU 
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: ZS 
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TD 4 
= .0 3 
R*SI 
T-I 
INT< 
RBS 



F 

R 

E CRMP 



I 

PI 



E>1 
TD 1 

INT 

S 1 

SI > 
INT 
TD 4 
INT 

fi" : 
SE. " 
INT 
PUT 

Z$: 

INT 



50 
5*B 
N < E 2 > 

X> 

< E X 1 THEN 180 

THEN 160 

40 

: : "TQTftL ROUNDS US 



1 8 T 
" N I C 
60 
"GD 
"REF 



HEN 450 

E SHODTINGl " 

BRCK TD FT. SI 
RESHER TRRININ 



"TRY fl 6 R I N ■; Y , N ) " 
TO BRS 



II y If 

"DK 



THEN 30 
RETURN 
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>RUN 
YDU ARE COMMRNDING fi GUN 
CREW. HIT WITHIN 100 YRRHS 
QF THE TARGET 

MAX RRMGE 41167 YDS 
TRRGET AT 17004 YARDS 
6UN ELEVRTIDN =77 
DVER BY 776 YARDS' 
GUN E L E V fl T I D N =78 
SHORT BY 534 YARDS 
GUN ELEVflTION =77.6 
TARGET DESTROYED! 
3 ROUNDS USED 

THE FDRUARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 21714 YARDS 
GUN ELEVATION =7 3 
DVER BY 1074 YARDS 
GUN ELEVATION =74 
SHORT BY 140 YARDS 
GUN ELEVATION =73.9 
TRRGET DESTROYED! 
3 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 16942 YARDS 
GUN ELEVATION =77.6 
TARGET DESTROYED! 
1 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 12531 YARDS 
GUN ELEVATION =82 
SHORT BY 1487 YARDS 
GUN ELEVATION =81 
SHORT BY 106 YARDS 
GUN ELEVATION =30.9 
TARGET DESTROYED! 
3 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 22631 YARDS 
GUN ELEVATION =73.1 
TARGET DESTROYED! 
1 ROUNDS USED 



TOTAL ROUNDS USED = 11 
NICE SHOOTING! 



TRY 

OK. 



AGAIN (.V ,tVj N 
RETURN TO BASE CAMP 



* * DONE * * 
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"I've just programmed our computer to give surprise birthday parties. " 

128 
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Potpourri 



Here are five programs that didn't seem to fit anywhere else. The first and 
the last are games, although you may come to think of the lunar landing 
simulation as a game also. One is a nifty simulation of smog, and the other 
calculates depreciation by three different methods. 
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Number Guessing Game 

Here is a computer program that plays the popular number guessing game. 
In it the computer picks a secret number between 1 and 100. You attempt to 
guess that number in as few tries as possible. 

There are many ways to go about guessing the secret number. Let some of 
your friends play this game and see what approaches they use to find the se- 
cret number. One approach is to start with a guess of 10. If this is too low, in- 
crease each guess by 10 until the computer says that a guess is too high. 
When this point is reached, start from the previous guess and increase each 
guess by 1 . This is the method used to solve some of the problems in the chap- 
ter on convergence. 

Another approach is to try to bracket the number between upper and lower 
limits and reduce the limits by steps until the number is finally found. Two of 
the convergence programs used this approach. 

Is one of these the best way? Well, these methods are not bad, particularly 
compared to starting with 1 and simply counting to 100 until the solution is 
found. But there is a better way. It is known as binary search. 

This technique involves dividing the search domain, in this case 1 to 100, in 
half, and then in half again, and so on until the secret number is found. Play 
the game many times using different approaches. In the long run you should 
find that the binary search approach is the most efficient. 

In Line 210, the program contains the statement that "you should not need 
more than 7 guesses." Why? If the secret number was between 1 and 128, 
what is the maximum number of guesses that would be necessary to find it? 
What if the number range were 1 to 130; then what would be the maximum 
number of guesses? 

Revise the program to choose a secret number between 1 and 10,000. Now 
the upper limit is 100 times the 1 to 100 game here which requires a maxi- 
mum of seven guesses to find the secret number; how many guesses will now 
be required? 

Can you write a program in which the roles of the computer and player are 
reversed? In other words, the computer will try to guess your secret number 
between 1 and some upper limit. After each guess, you enter L for low, H for 
high, or C for correct. Can you write this program so the computer can tell if 
you are cheating, i.e., giving it inconsistent clues? 
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IF N = 
IF G > 

PRI N 
6DTD 
PR IN 
eOTD 
PR IN 

ES" 
IF C 
PR IN 
GOTO 
IF C 
PR IN 
GOTO 
PR IN 

ED MO 



MIZE 

"I HRVE CHOSEN fl SE 
NTEGER BETWEEN 1 RND 
RV TD GUESS UHRT IT 
UILL GIVE YDU CLUES" 



hJT<:RND*100> + l 



T 

>3 

T 

2 
>7 

T 

*"> 

J- 

T 
RE 



YQ 
TH 
TH 
"T 

,0 
"T 

.0 
"C 

I T 
" Y 

■20 

T 

"G 

2 

"B 



UP. GUESS = 

EN 140 

E til 12 

D LOW" 

00 HIGH" 

DRRECT IN" 



HEN 180 
DU WERE 

HEN 210 
ODD JOB, 



C ; " G U 
LUCKY ! " 



UT YDU SHDULDH'T 
"THRN 7 GUESSES. 



PRINT " P L R Y R G fl I N < Y D R 



INPUT 
P R I tJ T 
IF flS = 
STOP 



flS 



THEN 20 



RUN 

I HRVE CHOSEN fl SECRET 
INTEGER BETWEEN 1 RND 100 
TRY TO GUESS WHRT IT IS. 
I UILL GIVE YDU CLUES 



YOUR GUESS- 
TDD L D IjJ 
'v' D U R GUESS 
TDD L D W 
Y D U R GUESS 
TDD L D W 
YDUR GUESS 
TDD LOW 
YDUR: GUESS 
CORRECT IN 
15 D n D J D B . 
P L R Y fl G R I N 

* * D D r J E * * 



= 50 

= 75 

= 87 

= 93 

= 97 

5 GUESSES 

< Y OR N ) ? 



131 



Depreciation — Three iVIethods 

This program shows how a piece of capital equipment depreciates accord- 
ing to three commonly used methods of depreciation: straight line, sum of the 
year digits, and double declining. 

The program asks for the original cost of the item, its expected life in years 
(the period of time over which it is to be depreciated), and its expected scrap 
(or sale) value at the end of that time. A table showing the annual depreci- 
ation for each of the three methods is then displayed. 



10 PRINT "PRDGRfiM CflLCULRTES 
DEPRE-" : "C IflTIDN BY THREE M 
E T H D D S " : " " 

2 INPUT " R 1 e I N ft L COST = " : C 
30 INPUT "LIFE DF ITEM CVEflR 
S ) = " : L 

40 INPUT "SCRAP VALUE =":S 
50 PRINT 

60 PRINT " STRAIGHT SUM DF 
DOUBLE" : "YEAR LINE DIGIT 
S DECLINING" s " " 
70 V=C-S 

8 D 1 = V / L 

9 Y = < C L + 1 ) .^ 2 > * L 
100 Z=L 

110 FDR X=l TO L 

12 D 2 = V * < Z ^' Y > 

130 Z==Z-1 

140 D3=2*C/L 

150 C=C-D3 

160 PRINT X5TABc:4>; 

170 Q=D1 

ISO GOSUB 2 80 

19 PRINT T fl B < 1 2 ) ; 

200 Q=D2 

210 GDSUB 280 

2 20 Q=D3 

230 PRINT TAB < 21)5 
240 GQSUB 280 

25 PRINT 

26 NEXT X 
270 STOP 

280 Q=IHTt:Q*100+.5>/100 

290 IF QMOOO THEN 310 

300 PRINT " "; 

310 IF Q>100 THE/«I 3 30 

32 PRINT " "? 

330 IF Q>10 THEN 350 

340 PRINT " " ! 

3 5 PRINT ; 
360 RETURN 
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PRDGRflM 
CIRTIQN 



CRLCULflTES DEPRE- 
BV THREE METHODS 



DRieiNflL CnST =8000 
LIFE DF ITEM <YEf=lRS> 
:SCRRP VALUE =500 



= 8 



STRAIGHT SUM DF 
YEAR LINE DIGITS 



1 
2 
3 
4 
5 
6 

7 
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9 37 
9 37 
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DDUBLE 
D E C L I N I H G 
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14 58 

12 5 

104 1 
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625 

4 16 

2 03 



67 

33 



6 

3: 

6 

3: 



2 
15 
1 1 25 
843 
632 
474 
3 55 
2 66 



( -J 
81 
61 
96 
97 



♦* DDME ** 



THIWkT 



=^ 




^ 




\((. 
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Smog Simulation 

This program is an adaptation of the smog model originally written by 
Herbert Peckham. The model assumes that vehicular traffic is the sole pro- 
ducer of smog, a somewhat poor assumption. It also assumes that most auto- 
mobile traffic occurs during the daylight hours and that traffic volume is very 
low (actually, zero) at night. The smog generated by the cars is dissipated by 
atmospheric conditions which vary depending upon sunlight, temperature, 
and weather. All of these conditions may be specified by the user. 

The model could be improved significantly by taking into account other 
sources of smog, by varying vehicular traffic according to the hour of the day, 
and by allowing daily variation of weather factors. Nevertheless, even in this 
rudimentary form it is interesting and instructive. 

A plot of the smog level is produced in Statement 370. Under some con- 
ditions, the smog level reaches a value that cannot be plotted because it is 
greater than the width of the screen (and printer). For runs with conditions 
like this, you might want to delete the plotting routine. A more elegant solu- 
tion would be to estimate the maximum value of the smog level from the in- 
put factors and calculate an appropriate plotting multiplication factor. 
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RINT 
RINT " 
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RINT 

RINT " D n 
RATE 



PRINT 
CLOUDY 

■I ■ 

INPUT 
PRINT 
PRINT 
OWN": 
PRINT 

cnoL 



TV CRRS PER RD 
S R N 6 E L E S 2 0" 

10 0": " T U L S fl 
: " S M D G CITY "5 



SMOG GENERRTID 
"19 5 R U T D .00 
D . I " : " B U S 
: " S M D G CITY " ? 



Y T I M E S M D 13 B R E R K 

IN PERCENT PER H 

CLERR .02" 

. 01 " s " SMDG 



R 1 

"N 
RR 



IGHTTIME SMDG BR 

T E f. V. PER H D U R > " 

HOT .05" 

.10":" S M D S 



I hJ P U T R 2 

PRINT 

PRINT "SMOG DISPERSION < 

R HDUR>" 

PRINT " H 

.00": "NO 

01 " : "SttQi 



I6H WIND RND RRI 
10 I N D AND DRY 
I G CITY 



190 

200 
EL 

210 

22 

230 

240 



250 

260 

270 

280 

290 

300 

310 

32 

33 


34 
350 
-.5 6 

37 

38 
390 
400 



INPUT R 3 
PRINT " 

PLOT 
R = R1 
K 1 = L 1 
T1 = INT<:C 
IF Tl/'2 = 

K 1=0 
R = R2 

S = S + K 1 * C 
IF S > ijr T 
T = T + 1 
T3=INT ( T 
T2=T-T3* 
J-i'INT OO 
IF T3.'-2 = 

PRINT T 2 
IF J< < 1 4 
X = 1 4 

PRINT T R 
GQTD 210 
PRINT T2 
GOTO 370 



HDUR Sl^nG LEV 



T - 6 ) -' 1 2 > 
INTCT1.^2)THEN 27 



*10-R*S-R3*S 
HEN 290 

/i2) 

12+1 

* £ + . 5 ■;■ -M 

INTtT3^2)THEN 39 

J " P M " 5 T R B ( 7 > 5 5<! 5 
THEN 370 

B < 1 3 + }< ~^ 5 " * " 

; " R M " ? T fl B ( 7 > ; X 5 
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Lunar Lander Simulation 

This program is one of the most popular computer simulations around. It 
is available in many versions; this one is adapted from the original program 
written in 1969. 

The program represents an exact simulation of an Apollo lunar landing 
module during the final descent. This portion of the descent would normally 
be controlled by the on-board computer backed up by another computer in 
the lunar orbiter, and still another computer on Earth. However, to exercise 
your knowledge of physics and to make an interesting game, we will assume 
that all three computers have had a simultaneous malfunction. Hence, it is up 
to you to land the spacecraft safely. 

To make a soft landing, you may change the burn rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (burn rate of 0) or of 
firing at a fuel rate of between 8 and 200 pounds per second. Engine ignition 
occurs at 8 pounds, hence values between 1 and 7 pounds are not possible. 
You have 16,500 pounds of fuel. This is 500 pounds more than an actual 
LEM has, which will give you a little margin for error. When you get pro- 
ficient, change Statement 80 to N = 16000 to simulate the real thing more 
closely. The capsule weight is 33,000 pounds. 

Not that this is the way to come in, but if you did not fire the rockets at all, 
the estimated time for a free fall descent is 120 seconds to impact (and a huge 
splat). 

Good luck! 
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Hammurabi 

Hammurabi is one of the all-time favorite computer games. On the one 
hand, it may be considered a game, but on the other it is an intriguing simula- 
tion of barter and management. 

Hammurabi is your servant as you try to manage the ancient city-state of 
Sumeria. The economy of the city-state revolves around just one thing — the 
annual crop of grain (probably soybeans). 

Each year, you must determine how many bushels of grain you wish to 
feed to your people (you'll quickly discover how much a person needs to sur- 
vive), how much you wish to use as seed in planting crops for the coming 
year, how much you wish to use for the purchase of additional land from 
your neighboring city-state, and how much you wish to put in storage. 

Of course, if you have a bad harvest or if rats overrun your grain storage 
bins, you may have to sell land in order to get enough grain to keep your peo- 
ple from starving, or to plant the land for the coming year. Unfortunately, 
disasters always seem to strike, forcing you to sell land, when the price is at 
an all-time low; but that's not any different from the real world. 

Most people start to play this game with noble ambitions. However, before 
long, they start longing for a plague to trim their growing population. Or they 
deliberately starve some people to keep things in balance (gosh, maybe these 
zero population growth people have something, after all!). 

Over the years, this game more than any other, has spawned a host of look- 
alikes, extensions, and modifications. Indeed, several manufacturers have 
taken my original with no changes whatsoever, put it in a fancy box, and 
charged a handsome price for it. Accept no imitations! Here is the original 
game (with the dialog shortened slightly) for you to run on your computer. 

If you want to experiment with changes, here are some suggestions. In the 
existing game, plagues randomly occur 15% of the time; lower this to 10% or 
5%. People now require a fixed amount of food; vary this amount slightly 
from year to year. Permit the construction of a rat-proof grain bin, but this 
must cost a fair amount. Introduce a mining industry as well as agriculture. 
How about fishing or tourism? Let your imagination run wild. Experiment! 
Have fun! 
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"Rats! A bacterium just ate the new micro-mini computer. 
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This Ideabook contains dozens of ways to make the most out of your 
computer for solving practical, everyday problems. The 50 ready-to-run 
programs demonstrate scores of different techniques for solving problems in 



PfKithematics, science and business. 



The ten chapters deal with solving problems by formulas and repetitive 



trials, convergence, recursion, compounding, probability, geometry, science. 



simulations, and drill and practice. 



Some of the problems demonstrote the capabilities of the computer; 
others identify its shortcomings. It is important to be familiar with both the 



strengths and weaknesses of your tools so you can recognize the types of 



jobs for which they are suitable. 
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